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Vorwort. 



Meine BemûhuDgen um die Verbesserung des mathematischen 
Apparats der experimentellen Psychologie (vgl. meine Studien 
zur Psychologie des Pessimismus, Wiesbaden 1904, mein Schopen- 
hauerbuch, Halle 1908, sowie die in den Sitzungsberichten der 
Wiener Akademie, math.-naturwiss. Klasse, 1915/16 verôfiFent- 
lichten Arbeiten „Bunte8te Reihen und Einge von Element- 
gruppen. Ein neues Problem der Kombinatorik", „Bunteste 
Ternen- und Qaaternenringe von harmonischer Struktur**, „Uber 
isonome harmonische Buntringe und eine merkwiirdige zwei- 
dimensionale Buntordnung'' usw.) ' fiihrten mich u. a. auch zu 
einer eingehenden Durchmusterung der Interpolationsmethoden. 

Dabei fiel mir auf, daB die psychologische Fachliteratur 
dièses Hilfsmittel recht stiefmutterlich behandelt. So bietet 
z. B. nur Alfred Lehmann in seinem vortrefiFIichen „Lehr- 
buch der psychologischen Methodik" (Leipzig 1906) eine aus- 
fiihrlichere Besprechung der Interpolation. 

Man stellt gewôhnlich die Forderung gleicher Intervalle. 
Demgegentiber scheint es mir wichtig,'[auf die Gaussische 
Interpolationsmethode hinzuweisen, deren Grundgedanke 
eine ôkonomische Ungleichteilung ist zum Zwecke môglichst 
vorteilhafter Berechnung von Mittelwerten. 

Nach einigen ersten Orientierungsversuchen bin ich zu der 
Uberzeugung gekommen, daB hier ein wertvolles Instrument zum 
Nachteil der Wissenschaft unbenutzt geblieben ist. 

Ich will nur ein Moment hervorheben, das sofort die hohe 
Bedeutung der Gaussischen Méthode fur psychologische Unter- 
suchungen ins Licht riickt. 

Jeder Experimentator muB darauf bedacht sein, bei der 
Versuchsperson nicht durch zu hàufige EingriflFe den natûrlichen 
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VI Vorwort 



Âblauf der BewuBtseinszustânde zu stôren. Da ist aiso eine 
Méthode, die mit einer geringeren Zabi von EiDgriffen dasselbe 
zu leisten verspricht, wie eine andere mit einer grôBeren Zabi, 
von vornberein als wertvoU zu begrûBen. 

Selbstverstândlicb beschrânkt sicb der Nutzen der Gaussi- 
scben Metbode nicbt auf dièses Anwendungsgebiet. Uberall, wo 
aus Sticbproben Durcbscbnittswerte zu ermitteln sind, kann sie 
gute Dienste leisten. 

Es erscbien mir zweckm'àfiig, zunâcbst die halbvergessene 
Gaussiscbe Abbandlung von 1814 mit ibren wicbtigsten Vor- 
arbeiten (Newton, Cotes) und der vereinfacbenden Darstellung 
Jacobis in deutscber Spracbe zugânglicb zu macben, damit 
das Interesse an diesem ausgezeicbneten matbematiscben Werk- 
zeug aucb auBerbalb des engeren Facbkreises neu belebt werde. 

Meinem Brader und matbematiscben Studienfreund Prof. 
Dr. 6. Kowalewski in Prag danke icb fur die Mitwirkung an 
der Korrektur und ftir verscbiedene facbmânniscbe Auskûnfte. 

Aucb môcbte icb nicbt unterlassen, das bocbberzige Ent- 
gegenkommen, das mir der Herr Verleger durcb Ubernabme 
dieser Scbrift in einer ftir den Bucbbandel so tiberaus scbweren 
Zeit bewiesen bat, riibmend bervorzubeben. 

Kônigsberg i. Pr., im September 1916. 

Arnold Kowalewski. 
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Die Differentialmethode. 

(Von Isaac Newton, 1711.) 



Brste Anlstellnng. 

Die Âbszisse einer krummlinigen Figur setze sich 
aus irgendeiner gegebenen Grôfie A und einer un- 
bestimmten GrôSe x zusammen, die Ordinale bestehe 
aus irgendeiner Ânzahl gegebener GrôSen b, e, d, e..., 
die mit ebensovielen Oliedern der geometrischen Pro- 
gression X, x'^, x\ x\ . . . der Reihe nach multipliziert 
sind, und zu gegebenen Âbszissenpunkten errichte man 
die entsprechenden Ordinaten. Dann sage ich, daB 
man die ersten Differenzen der Ordinaten durch ihre 
Intervalle dividieren kann, die Differenzen der so di- 
vidierten Differenzen durch Paare von Ordinateninter- 
vallen, die Differenzen dieser so dividierten Diffe- 
renzen durch Ternen von Ordinatenintervallen und so 
fort ins Unendliche. 

Man setze in der Tat fiir den unbestimmten Teil x der 
Abszisse nacheinander beliebig gegebene GrôSen p, q^ r^ s, i,, . . 
und errichte in den Endpunkten der so gegebenen Abszissen 
die Ordinaten a, fi, y, S, b, ... Dièse Abszissen und Ordinaten, 
die Ordinatendififerenzen dividiert durch die Abszissendifferenzen 
(die eben die Ordinatenintervalle sind)^ die Differenzen dieser 
Quotienten dividiert durch je zwei aufeinanderfolgende Ab- 
szissendifferenzen * usw. sind in folgender TabeUe verzeichnet**: 



* Im Original steht irrtUmlich ,,Ordinatendifferenzen'^ 
** Im Original lautet das konstante Glied der Ordinate J» Wir 
haben es in a verwandelt, weil schon das konstante Glied der Abszisse 
mit A bezeichnet wird. Beide gleichzusetzen, ist Willkiir, 
Kowalewskl, Newton oBw. 1 
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AbszisBen 



Ordinaten 



A + p 


a + 6^ + Gp* + dp^ + ep* = a 




A-h q 


a + 6 (/ + c fy* + rf ç' + e î* = /? 




A +r 


a + ôr + cr* + rfr*+er* = y 




A ■{- s 


a + ôs + es* + rfs' + e«* « (^ 




A + t 


a + 6'^ + c^« + d^» + e^* = e 




Reihe der 
divid. DifiF. 


Divisionsergebnisse 





« — 


^ 


P- 


9 


ti- 


r 


9- 


r 


r - 


ô 


r — 


s 


ë - 


e 



s - t 



b + c ip + q) ^ d {p^ -^ p q -^ q^) + e {p^ ■{- p^ q + pq* + q*) = C 
b -{■ c {q + r) + d{q^ -h q r + r*) + eiq^ + q^r + q r^ + r^) = 7f 
b + c(r + s) + d{:r^ + r s + s^) + e {r^ -{- r* s + rif* + s^) = & 
b + cis + t) + d {s* + s t + i^) + e {s^ -{- sU + s i^ -h t^) ^ x 



y> - r 

q-s 
it - X 

r - t 



c -{- d{p + q + r) -h e{p* -\- p q + q^ -i- p r + qr + r^) =^ X 
c + d{q + r -{■ s) + e{q^ + qr 4-r* + g'S + rs + 5*; = /ii 
c -^ d{r + s + t) -h e(r* +rs'h8'-^rt+st'\-t*)^v 



p — s 

A^ - y 

q-t 



d-\-e(p-hq + r-hs) = è 
rf + c(î + r + s + t) ^ n 






Zweite Ânfstellnng. 

Uûter denselbea Voraussetzungen und bei einer 
endlichen Anzalil von Gliedern b, c, d, e, ... behaupte 
ich, daB der letzte Quotient gleich dem letsiten der 
Glieder b, c, d, e . . . sein wird, daB ferner dùrch die 
ubrigen Quotienten^ die iibrigen Glieder b, c, d, e ..• 
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bestimmt sein werden und mit diesen die krumme 
Linie parabolischer Natur, die durch die Endpunkte 
aller Ordinaten hindurchgeht. 

In der Tat war in der obigen Tabelle der letzte Quotient a 
gleich dem letzten Gliede e. Multipliziert man dièses Glied mit 
der gegebenen Summe p + q + r + s und subtrahiert es von 
dem Quotienten |, so bleibt das vorletzte Glied d ûbrig. Nimmt 
man femer die bereits bekannten GrôSen 

dij> + q + r) + e{p^ + pq + q^ +pr + qr + r^) 

von dem Quotienten k fort, so bleibt das Yorvorletzte Glied o 
tibrig. Nimmt man weiter die bereits bekannten Grôfien 

c{p + q) + d{p- + pq + q*) + e{p^ + p^ q + pq^ + q^) 

von dem Quotienten Ç fort, so bleibt das Glied b, Durch âhn- 
liche Bechnung wûrde man, wenn mehr Glieder da wàren, aile 
durch ebenso viele Ordnungen von Quotienten gewinnen. Werden 
dann die gegebenen Grôfien 

bp + cp^ + dp^ + ep^ 

von der ersten Ordinate a abgezogen, so bleibt das erste Glied a 
der Ordinate* tibrig, und die GrôBe** 

a + bx + cx^ + dx^ + ea;* 

steUt die Ordinate der Kurve parabolischer Natur dar, die durch 
die Endpunkte aller gegebenen Ordinaten hindurchgeht, wenn 
^ + a; die Abszisse ist. 

Ans diesen Aufstellungen lâfit sich das Folgende leicht 
herleiten. 

Dritte Anfstellung. 

Man telle eine beliebige Gerade*** ^^^ in beliebig 
viele gleiche Teile A^A^, A^A^, ^8-^4» -^^A» • • • ^^^ ®^" 
richte in den Teilpunkten die Parallelen A^B^, A^B^, 
A^B^y ... Es soll eine geometrische Kurve parabolischer 
Natur gefunden werden, die durch die Endpunkte B^, 
i^2, jBg, ...aller errichteten Strecken hindurchgeht. 

* Das Original hat: „A der Abszisse". 
•* Das Original hat: ,,4 + ôa; + ex* + . . ." 
*** Newton lâfit den Index 1 immer fort. 

1* 
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Man sache von den errichteten Strecken A^B^, A^B^, 
^3^3, .. . die ersten Diflferenzen \^ b^, èj, . . ., die zweiten c^, 
Og, C3, . . ., die dritten d^, d^, rf,, ... usf., bis man zur letzten 
Differenz gelangt, die hier i^ ist. 

Dann nebme man, mit der letzten Differenz beginnend^ in 
den Eolonnen oder Ordnungen der Differenzen abwécbselnd die 



9, 



f 



^. 



fs 



^. 




Fig. 1. 

mittelsten Differenzen und die arithmetischen Mittel ans den 
beiden mittelsten. So fahre man ordnungsmâfiig fort bis zur 
Reihe der ersten Glieder -Ij-B^, A^B^^ A^B^, ... Auf dièse 
Weise ergebe sich: 

k, Ij m, w, 0, Py q, Tj s, . . .^ 

\i'0 das letzte Glied die letzte Differenz bezeichnen môge, das 
Yorletzte das arithmetische Mittel ans den beiden vorletzten 
Differenzen, das vorvorletzte die mittelste der drei vorvorletzten 
Differenzen usw. bis zum ersten, welches das mittelste der Glie- 
der* A^B^, A^B^j A^B^, ... sein wird oder das arithmetische 
Mittel ans den beiden mittelsten. Das erstere ist der Fall bei 



Das Original hat irrtûmlich: ^^A^ A^^ A^". 
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uDgerader, das zweite bei gerader Anzahl der Glieder A^B^, 

Erster Fall. 
Im ersten Falle sei A^B^ jenes mittelste Olied, also 

Es werde nun die Ordinale PQ errichtet, und man setze 

Fa£t man dann die Glieder folgender Progression 

X X x^ -^ 1 x_ a;* — 4 x a;^ — 9 



iC 



a;* - 16 a; a;* - 25 a; a;* - 36 



T' 9x ' 10' lia; ' 12' 13a; ' * ' ' 

fortlaufend zu Produkten zusammen, so entstehen die Glieder 

- x* x^ — X X* — x^ a;* — 5 a;* + 4 a? 

^' ^' ~2"' 6 ' 24 * 120 ' 

a;0 - 5 a;* + 4 a;* a?^ - Ux^ -i- 49a;» - 36a? 

720 ' 5040 ' * • • 

Wenn man mit ihnen die bezUglichen Glieder der Reihe 
k, ly m, n, Oj p, , . . 
multipliziert, so wird das Aggregat dieser Produkte 

, , , , a;^ , x^ — X , X* — X* , a;* — 5 a?» + 4 a: 

k + xl+—m+—^—n + ^^^o + — P+ ... 

die Lange der Ordinale PQ sein.^ 

Zweiter Fall. 

Im zweiten Falle seien A^B^, A^B^ die beiden mittelsten 
Glieder, also 

2 ^' * ^' — 2 — "^ ' 
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«3 — *•> 2 "" ^ > h ^ Pj 2 ~ ^' ^1 '^ • 

Es werde nun die Ordinate PQ errichtet und man halbiere 
A^A^ in O und seize OP=^x, FaBt man dann die Glieder 
folgender Progression 

^«1 t ^ , 25 . 

i ^ 4^ x_ 4^ jç^ 4 

'1' 2a; ' 3 ^ 4a; ' T' 6« ' 

a;« ;,- 

X 4 



7 ' 8a; ' ••• 

fortlaufend zu Produkten zusammen, so entstehen die Glieder 

1 4a;« - 1 4a;^-a; 16 a;* - 40 a;» + 9 

' ^' 8 ' 24 ' 384 '••• 

Wenn man mit ihnen die beziiglichen Glieder der Reihe 
k, l, w, w, 0, p, q, . . . 

multipliziert, so wird das Aggregat dieser Produkte 

, , , , 4a;*-l , 4a;»-a; , 16a;* - 40a;* + 9 

k + xl + —^—m+^^^n + 3gj 0+ ... 

die Lange der Ordinate PQ sein.^ 

Hierbei ist aber zu bemerken, dafi die Intervalle A^A^, 
A^A^, A^A^, ... als Einheiten vorausgesetzt werden, und dafi 
die Differenzen durch Abziehen der unteren Grôfien von den 
oberen gebildet werden miissen. • Man hat also A^B^ von A^B^, 
A^B^ von A^B^, b^ von b^ usw. abzuziehen und A^B^^ A^B^^^b^, 
A^B^ — -4g 5^3 = ^2, ^1 — ^2 ~ ^1 ^^^' ^^ setzen. Uberdies mufi 
man, wenn sich auf dièse Weise négative Differenzen ergeben, 
die Zeichen ândem. 

Vierte Anfstellnng. 

Man teile eine Gerade in beliebig viole ungleiche 
Telle A^A^, ^^3* ^z^a^ ^aA> • • • ^^^ errichte in den 
Teilpunkten die Parallelen A^B^, A^B^, -^gfîg, ... Es 
soll eine geometrische Kurve parabolischer Natur ge- 
funden werden, die durch die Endpunkte B^, B^y B^, ... 
aller errichteten Strecken hindurchgeht.* 



* Im Original steht irrtumlich: e^ + e^ ohne den Nenner 2. 
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Es seien die Punkte ^,, B^, 2?g, B^, B^, B^\ B^, ... ge- 
geben, und man fàlle senkrecht auf irgendeine Abszisseaachse 
A^Aj die Ordinaten B^A^, ^2^2» • • • 









1: 




]': 








A"' 






9, 


l' 


9, 


7 


7<^ 


'O, 



Fig. 2. 



Man seize 
Al Bi — A2 B^ ^ , A ^B^ — Af^B^ __ - 



Ai B4 — A^ B» 



= \, 



/l,B, - M.B, 



darauf 



= ^. 



A ^B^^ AjB, 
AqA^ 






*' - *» = c 



6. -ft, _ 



weiter 



und 



'~^d„ 



A, A, 



^ - <•« 



= «^2» 



«1^ 



AfAg 



= e. 



2> 



A, A, 

rfg — d^ 



^3» 



^ = ^3, 



^3^7 



««. 



So fahre man fort bis zur letzten DiflFerenz. 

Nachdem man in dieser Weise die DiflFefenien gebildet und 
aie durch die Ordinatenintervalle dividiert bat, nehme man in 
ibren Kolonnen oder Reiben oder Ordnungen von der letzten 
beginnend abwechselnd die mittelsten und die aritbmetiscben 
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Mittel aus den beiden mittelsten und fahre fort bis zur Beihe 
der ersten Glieder A^B^, A^B^, ... Auf dièse Weise >er- 
gebe sich 

k, l, m, n, 0, p, q, r, . . ., 

wo das letzte Glied die letzte Dififerenz bezeichnet, das vor- 
letzte das arithmetische Mittel aus den beiden vorletzten, das 
vorvorletzte das mittelste der drei vorvorletzten usw. Das erste 
Glied k wird die mittelste Ordinate sein, wenn die Anzahl der 
gegebenen Punkte ungerade, oder das arithmetische Mittel aus 
den beiden mittelsten, wenn ihre Anzahl gerade ist 

Erster Fall. 
Im ersten Falle sei A^B^ die mittelste Ordinate, also 

Es werde die Ordinate PQ errichtet und auf der Basis A^A^ 
ein beliebiger Punkt angenommen. Man setze 

OP^x 

und fasse die Glieder folgender Progression 

OAs + OA^ (X - OA^) (X - OA^) 



1, aj — OA^j 



2 ' x-\(OA^-¥ OA^) 



_ OA^ + OA^ Og-Q^,)(a;-QJe) _ O^i + OA, 

^ 2 ' X'-\(OAt+ OA^y ^ 2 

fortlaufend zu Produkten zusammen und halte die s6 eût- 
stehende Progression fest; oder, was dasselbe ist, man fasse, die 
Glieder folgender Progression 

1, x^ OA,, {x - OA,){x - OA,), {x - OA^){x - OA,) 

iortlaufend zu Produkten zusammen und multipliziere die so 
entstandenen Glieder mit den entsprechenden Gliedern fol< 
gender Progression: 

. ^ OA, + OA, ^ OA, 4- OA^ ^ OA, + OA, 

A, X - , X - , X— 2 
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Newton: Die Differentialmethode 9 

Dann entsteben die unpaarigen Glieder, wSlhrend die yoUstân- 
dige Progression lautet 

{x^OA^){x^OA^){x^OA^), ... 

Oder man seize 

OA^ = a, OA^^p, OA,^ry OA^^S, OA,^ej 

OA^ = f, OAj^ 17; 

OA^ + OA ^ 0A^+ OA^ ^ OA^ + O Ay . 
-—"2 =^> 2 ^^' 2 ^ 

Und aus der Progression 

1, x — S, (x — y)(x -^ b), (« - /S)(a; — Ç), 

sammle man die Glieder, die mit 

1, 05 — i9-, X -- Xf X " X 

multiplizieit werden. Dann fUge man die fehlenden Zwischen- 
glieder ein, wobei sich die yoUstândige Seihe 

1, x-d, x^^{3+&)x + d&, 
x^'-'{3 + 2&)x^ + {ye + 2S&)x--y3e, . .. 

ergibt. Mit ihren Gliedem ist die Beihe k, l, m, n, 0, . , . zu 
multiplizieren. Das Aggregat der Produkte 

k + [x^3)l + [x^ -{3 + &)x + 3&']m + ... 

wird die Lange der Ordinate PQ sein. 

Zweiter FalL 

Im zweiten Falle seien A^B^^ ^s^t ^^ beiden mittelsten 
Ordinaten, also 

— . 2 = «, 0^ = i, — - — = m, rfj = n, 

— 2 — ' /2 =-P> 2 — "^^^ ^3=r. 

Die Koeffizienton der unpaarigen k, m, 0, q, . . . werden 
durcb fortlaufende Multiplikation der Glieder folgender Pro- 
gression entsteben: 

1, {x - OA^){x - 0^5), {x - OA^){x - O^g), (x - OA^){x - O^^). 
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10 Newton: Die Differentialmethode 

Die Koeffizienten der tibrigen ergeben sich durch Multipli- 
kation dieser mit den Gliedem folgender Progression: 

OA^-^-OA^ OA^JrOAf, OA^Jr OAi 

^ 2 ' 2 ' 2 ' 

2 

d. h. es wird 

die Ordinale PQ sein, oder 

PQ = k'\'X \l +x \ * X \m + x \ • X \ * x \n+... 

Oder man setze 

r (a; - ^^1±^M = v, r(rr - OA,)(a; - OA,) = (p, 

^r — ^^ — ^j"^- 

Dann wird 

sein. 

Fûuf te ' Au! stellung. 

Wenn irgendwievieleGlieder einer beliebigenReihe, 
auf bestimmte Intervalle verteilt, gegeben sind, ein 
beliebiges Zwischenglied môglichst angenàhert zu 
finden. 

Auf einer Geraden von gegebener Lage errichte man 
senkrecht die gegebenen Glieder, getrennt durch die gegebenen 
Intervalle. Durch die Endpunkte der Lote werde nach den 
vorigen Aufstellungen eine krumme Linie parabolischer Natur 
gelegt. Sie wird aile Zwischenglieder die ganze Reihe hin- 
durch treffen. 
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Sechste AnfsteUang. 

Eine krummlinige Figur, von der einige Ordinaten 
gefunden werden kônnen, môglichst genau zu qua- 
drieren. 

Man ziehe durch die Endpunkte der Ordinaten eine Kurve 
parabolischer Natur mit Hilfe der vorhergehenden Aufstellungen. 
Sie wird nâmlich eine Figur begrenzen, die immer quadriert 
werden kann, und deren Flâcfae ziemlich genau gleich der 
Flâche der vorgelegten Figur sein wird. 

Anmerkung. 

Dièse Aufstellungen sind nûtzlich zur Herstellung von Ta- 
feln durch Interpolation von îleihen, sowie zur Lôsung von 
Problemen, die von der Quadratur der Kurven abhàngen, be- 
Bonders wenn die Intervalle der Ordinaten klein und unterein- 
ander gleich sind und man fiir jede gegebene Anzahl von 
Ordinaten Regeln herstellt und fur den Gebrauch aufbewahrt. 
Liegen z. B. vier Ordinaten in gleichen Intervallen vor, so sei 
A die Summe aus der ersten und vierten, B die Summe der 
zweiten und dritten und R das Intervall zwischen der ersten 
und vierten. Dann wird die neue Ordinate in der Mitte aller 
9J9- A 
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sein und die ganze Flâche zwischen der ersten und 
vierten • R. 

o 

Man beachte noch folgendes. Wenn die Ordinaten gleiche 
Abstânde voneinander haben und man die Summen der Ordi- 
naten nimmt, die von der mittleren Ordinate gleich weit ent- 
fernt sind, und das Doppelte der mittleren Ordinate, so wird 
eine neue Kurve zusammengesetzt, deren Flâche durch weniger 
Ordinaten bestimmt wird und gleich der Flâche der. vorigen 
Kurve ist, die man finden soUte. Ja man kann auch als 
neue Ordinaten die Summe der ersten und zweiten Ordioate, 
die Summe der dritten und vierten, die Summe der fûnften 
und sechsten nehmen usw., oder die Summe der drei ersten 
Ordinaten, die Summe der drei nâchsten, die Summe der drei 
weiter folgendeu, oder die Summe von je vier oder die von je 
fûnf Ordinaten. Die Flâche der neuen Kurve wird gleich sein 
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der Flâche der zuerst vorgelegten Kurre. Und so wird, wenn 
man irgendeine Anzahl von Ordinaten der zu quadrierenden 
Karve bat, die Qaadratur derselben auf die Quadratur einer 
anderen Kurve durch weniger Ordinaten zurûckgefiihrt werden. 
Darch irgendeine Anzahl gegebener Punkte lassen sich aber 
nicht nur Eurven von paraboliscber Natur, sondern auch un- 
zâblige andere von verscbiedenen Arien legen. 

Es seien CDE, FOH zwei Kurven mit der gemeinsamen 
Abszisse AB und den aaf derselben Geraden liegenden Ordi- 
naten BD, BQj und die Beziehung 
zwischen diesen Ordinaten werde 
durch irgendeine Gleichung definiert. 
Es môgen Punkte in beliebiger An- 
zahl gegeben sein, durch welche 
die Kurve CDE hindulrchgehen muB. 
Durch jene Gleichung werden dann 
ebenso viele neue Punkte sich er- 
Fig. 3. geben, durch welche die Kurve 

FOH hindurchgehen wird. Man 
beschreibe mittels der friihere'n Aufstellungen die Kurve FOH 
von paraboliscber Natur, welche durch aile jene neuen 
Punkte geht. Dann wird durch dieselbe Gleichung die Kurve 
CDE geliefert werden, die durch aile zuerst gegebenen Punkte 
hindurchgeht 




Uber die Newtonsche Differentialmethode. 

(Von Roger Cotes, 1722.) 



Unser beriihmter Newton lehrt gegen Ende des dritten 
Bûches seiner Prinzipien der Naturphilosophie, wie man 
aus irgendeiner Anzahl beobachteter Orter eines Kometen den 
Ort desselben fur irgendeine Zwischenzeit findet. Bei dieser 
Gelegenheit beschreibt er jene berûhmte und sehr allgemeine 
Méthode, durch beliebig viele gegebene Punkte eine 
geometrische Kurve zu ziehen^ eine Méthode, die er ein- 
mal in einem Briefe an Leibniz mit Becht zu den schônsten 
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Errangenscbaften z&hlte.^ Dièse werde ich nun, da sie von 
gro£er Wichtigkeit besonders in astronomischen Dingen ist und, 
soviel ich weifi, bis jetzt noch Biemand eine Begrûndung des 
Yerfabrens entwickelt bat, im folgenden etwas ausfûbrlicber dar- 
iegen und nocb einiges von mir selbst binzufûgen. 

Ërste Anfstellnng. Lehrsatz. 

Es sei eine beliebige Reibe von GrôBen vorgelegt 
Wenn man aile Produkte von irgendeiner Dimension, 
die sich ans allen, mit Ausnabme der letzten, bilden 
lassen, zu einer Summe vereinigt und ebenso aile Pro- 
dukte derselben Dimension, die sicb aus allen, mit 
Ausnabme der ersten^ bilden lassen» zu einer Summe 
vereinigt, um dann die zweite Summe von der ersten 
abzuzieben und den Best durcb die Differenz zu divi- 
dieren, die durcb Subtraktion der letzten der vor- 
gelegten Grôfien von der ersten entstebt, so bebaupte 
icb, daB die durcb solcbe Division sicb ergebende 
GrôBe gleicb ist der Summe aller Produkte von der 
nâcbst niedrigeren Dimension, die sicb aus der ganzen 
Keibe der vorgelegten GrôBen bilden lassen. 

1. Beispiel Wenn die vorgelegte Reibe nur aus zwei 
GrôBen a und b bestebt, und die aus diesen zu bildenden Pro- 
dukte von vierter Dimension sind, so ist das einzige Produkt 
aus allen^ mit Ausnabme der letzten, d. h. aus der ersten, a\ 
und das einzige Produkt aus allen, mit Ausnabme der ersten, 
d. b. aus der letzten, b*. Zieht man dies von dem ersteren ab, 
80 wird sicb der Rest a* — ô* ergeben. Dividiert man ibn 
darcb a — b, so erbâlt man 

a^ + a^b + ab^ + b^, 

also wirklicb die Summe aus allen Produkten dritter Dimen- 
sion, die sicb aus der ganzen Reibe der GrôBen a und b bilden 
lassen. 

2, Beispiel. Die vorgelegte Reibe bestebe nunmebr aus 
drei GrôBen a, b, o und es seien die aus ibnen zu bildenden 
Produkte von dritter Dimension. Die Summe der Produkte aus 
allen, mit Ausnabme der letzten, wird also sein 

a^ + a^b + ab^ + b^ 
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und die Summe der Produkte aus allen, mit Ausoahme der 
ersten, wird sein 

h^ + b^c + bc^ + cK 

Man subtrahiere dièse Summe von der ersteren. Dann bleibt 
der Best 

flS ^ ctH + ab^ - b^c - bc^ - c^. 

Dividiert man ihn durch a — c, so erhâlt man 

a^ + ab + ao + b^ + bc + c^, 

d. h. die Summe aller Produkte zweiter Dimension, die sich 
aus der ganzen Reihe der GrôBen a, b, o bilden lassen. 

3. Beispiel. Es werde vorgelegt eine Eeihe von vier GrôBen 
a, 6, c, d, und die aus diesen zu bildenden Produkte seien von 
zweiter Dimension. Die Summe der Produkte aus allen, mit 
Ausnahme der letzten, wird also sein 

a^ + ab + ac + b^ + bc + c^, 

die Summe der Produkte aus allen, mit Ausnahme der ersten, 
dagégen 

h^^hc+bd-\-c^ + cd + dK 

Man subtrahiere dièse- Summe von der ersteren. Dann bleibt 
der Rest 

a^ + ab + ac --bd — cd ^ d^. 

Dividiert man ihn durch a — d, so erhâlt man 

a + b + c + d, 

also die Summe aller Produkte erster Dimension (wenn man so 
sagen darf), die sich aus der ganzen Reihe bilden lassen. 

4. BeispieL Wenn eine Reihe von fûnf GrôBen a, % c, rf, e 
vorgelegt ist und die aus ihnen zu bildenden Produkte von 
erster Dimension sind, so wird die Summe der Produkte aus 
allen, mit Ausnahme der letzten, 

a + b + c + d 

sein und die Summe der Produkte aus allen, mit Ausnahme 
der ersten, 

b + -h d + e. 

Subtrahiert man die zweite von der ersten, so bleibt der Rest 
a — e. Dividiert man ihn durch a — e, so erhâlt man 1, das 
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eiozige Produkt nullter Dimensioo, das sich aus der ganzen 
Eeihe bilden lâBt. 

Zweite Aufstellnng. Lehrsatz. 

Die unbestimmte Gleichung 

y = f+gx + hx^ + kx^ + . . ., 
von der aile Koeffizienten fj g, h^ k, . . . gegeben sind, 
werde so begrenzt, daB sie sicb auf irgendeine gege- 
bene analytische Ordnung bezieht. Fiir x setze man 
der Reihe nach irgendwelche verschiedenen Werte a, b, 
c, d, e, , , , ein, und es môgen sich dabei ftir ^ beziiglich 
die Werte a, /9, ;', Ô, «, . . . ergeben. Man stelle dièse 
zusammengehôrigen Werte in einer Tabelle einander 
gegeniiber: 

'la « ' 



U8W. 



und fûge die Differenzen hinzu, auf solche Weise, daB 

// = a — b, X" = b — c, i!'' = c — d^ X"" =i d — e, ... 

fi =s a — G, fi' =^ b -- d, fx" = c — e, . . ^. 
V =^ a — d, 1/" = è — 6, . . . 

ist usw. Ferner sei 

r'^a^^, r" ^§^r, r'" ^ y ^ 8, r""=:^-€, . 

und dann 

Wiederum sei 
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und dann 






So fahre man fort. 
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Dann werden die so diyidierten ersten Differenzen 
Qi q'\ i/"j • • • einander und dem Koeffizienten hôchsten 
Grades g gleich sein, wenn die Gleichung von erster 
Ordnung ist, y^f + gx. Es werden die zweiten g\ a% 
a"\ ... einander und dem Koeffizienten hôchsten Gra- 
des h gleich sein, wenn die Gleichung von zweiter Ord- 
nung ist, y = /' + ^a5 + Aa;*. Die dritten t', t", t'", . . . 
werden einander und dem Koeffizienten hôchsten Gra- 
des h gleich sein, wenn die Gleichung von dritter Ord- 
nung ist, y ^f -^ gx'\- hx^ •\''koi?. So kann man fort- 
fahren, von welcher Ordnung auch die Gleichung 
sein mag. 

1, FalL Es môgen in der vorgelegten Gleichung aile Glie- 
der fehlen, his auf das vom hôchsten Grade, dessen Koeffizient 
der Einheit gleich sei. Dann wird y gleich einer reinen Potenz 
von X. Zq welcher Dimension in a, b, c, d, e, . . ., daher auch 
a, /S, ;', S, e, . . . aufsteigen môgen, so ist (nach der 1. Aufstel- 
lung) klar, daB ()' die Summe aller Produkte der nâchst niedri- 
geren Dimension sein wird, die sich aus den Grôfien a und b 
hilden lassen, ç" die Summe aller ans b und c, ç'" die Summe 
aller aus c und d, q'" die Summe aller aus d und e. Âhnlich 
wird a' die Summe aller Produkte der nachst niedrigeren Dimen- 
sion sein, die sich aus den drei GrôBen a, 5, c bilden lassen, 
a" die Summe aller aus 6, c, d, ferner a" die Summe aller aus 
e , d, e. Wiederum werden t', t" die Summen aller Produkte 
noch niedrigerer Dimension aus je vier GrôBen sein. So geht 
es weiter, bis mad schlieBlich auf Einheiten gekommen ist, wie 
bei dem vierten Beispiel der ersten Aufstellung. Dies muB ent- 
weder bei den ersten dividierten Differenzen eintreten oder bei 
den zweiten oder bei den dritten usw., je nachdem y gleich x 
ist oder gleich dem Quadrat oder dem Kubus von x usw. 

2. FaU. Unter denselben Voraussetzungen sei nun der 
Wert jenes hôchsten Koeffizienten irgendein anderer. Dann 
wird y gleich einer Potenz von x sein, multipliziert mit jenem 
Werte. Es werden somit a, /?, y, S^ e, . . . bezûglich gleich 
derselben Potenz von a, b, c, d, e, . , . , multipliziert mit dem- 
selben Koeffizienten, und jener Koeffizient wird den ersten und 
den folgenden Differenzen ' anhaften. Daher werden jene letzten 
Differenzen, die im vorigen Falle einzein der Einheit gleich 
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waren, im vorliegenden Falle einzeln diesem Eoeftizienten 
gleich sein. 

3. FaU. Es môgen jetzt beliebig viele Glieder in der Glei- 
chung vorhanden sein. Ofifenbar werden jene Bestandteile der 
Differenzen zaerst gleich werden und verschwinden, die ans den 
niedrigeren Gliedern, dann jene, die aus den hôheren hervor- 
gehen, nnd zuletzt bleiben nur jene tibrig, die dem hôchsten 
Gliede entstammen. Es werden daher jene letzten Difierenzen 
auch in diesem Falle einzeln dem Koeffîzienten des hôchsten 
Gliedes gleich sein, was zu beweisen war. 



Dritte Anfstellnng. Problem. 

Man denke sich unter denselben Voraussetzungen 
wenigstens so yiele Werte, a, 6, c, d, ... der unbestimmten 
Grô£e x gegeben, als die vorgelegte Gleichungsordnung 
Glieder baben kann. Es seien ferner ebenso YÎele a, /?, 
^, J, . . . gegeben als zugehôrige Werte det unbestimmten 
GrôBe y. Wâhrend nun die Gleichung noch unbekannt 
ist, sei verlangt, den Wert J von y zu finden, wenn x 
als Wert von x angenommen wird. 

Man fiihre ailes nach dem Muster der vorigen Aufstellung 
aus, bis man schlieBlich erreicht, dafi gewisse dividierte Diffe- 
renzen, etwa t\ t", . . . gleich sind. Dann nehme man X = 
«— o, fi ^ z ^ b, V = z -^ e usw. 
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Man mufi so weit gehen, bis ebenso viele X, fi, v, . 
Sind, als q\ a, t', ... Dann wird sein * 

was gesucht war. 

kowalewski, Newton usw. ^^ 
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Man bedenke in der Tat, daB 

Dann werden (nach der 2. Aufetellung) t, t', t" aile einander 
gleich sein. Daher wird sein 

t = Tp =z T p nnd a ^ a' + t = a + TP 
nnd 

s =^ a fi = GfA + T fip und ç ss q' ^ s ^ ç' + afi + r'/i v 

nnd 

endlich 

was za beweisen war. 

Yierte Anlstellnng. Problem. 

Bei denselben Voraussetzungen und Daten soll die 
unbekannte Gleichung ermittelt werden, und zwar 
durch Bestimmung der Koeffizienten f, g, h, k, . . . der 
einzelnen Glieder. 

Man seize z = 0. Dann wird Ç (welches nach der yorigen 
Aufstellung zu finden ist) gleich der GrôBe f sein; wenn nâm- 
lich z oder x verschwindet, so verschwinden zugleîch die Glie- 
der gx + hx^ + kx^ + ... und es bleibt f gleich y oder Ç. 
Nun ziehe man die Bekannte f von beiden Seiten der Glei- 
chung ab und dividiere, was Ubrig bleibt, durch x. Es ergibt 
sich 80 

^'' ^ =^g + hx + kx^+ ... 
Man seize 

^ X 

Dann bat man die niedrigere Gleichung 

y=g + hx + kx^+ ... 

Ebenso schreibe man fur a, /S, y, S, s, . . . in der Tabelle ce' 
/5', /, ^, «',..., so daB 

a ^— , fF - j , . .. 
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ist. Nun setze man von neuem x ^ 0. Dann wird f (welches 
nach der vorigen Aufstellung zu finden ist] gleich der Grôfis g 
sein ans den oben angegebenen Griinden. Durch Wiederholung 
dièses Verfahrens wird man die Gleichung fortgesetzt herab- 
drûcken, bis schliefilich aile Eoeffizienten bestimmt sind, wie 
es verlangt war. 

Dasselbe auf andere Weise. 

Es steht nach der 2. Aufstellung fest, daB die letzten divi- 
dierten Differenzen gleich dem Eoeffizienten k des hôchsten 
Gliedes sind. Daher wird jener Koeffizient k bekannt werden, 
indem man nach der oben beschriebenen Méthode die Diffe- 
renzen nimmt. Nun ziehe man das hôchste Glied kx^ von 
beiden Seiten der Gleichung ab. Dann erhâlt man 

y -- kx^= f+ gx + hxK 

Man setze y = y — kx^. Dann bat man die niedrigere Gleichung 

y' ^ f+gx + hxK 

Ebenso schreibe man fur a, jS, /, 5, . . . in der Tabelle a\ /9', 
/, ^, . . . , 80 daB a = a — ka^, ^ = ^ - kh^, . . . ist. Man 
nehme wiederum die Differenzen, um den Eoeffizienten h des 
hôchsten Gliedes dieser néftien Gleichung zu erhalten. Durch 
Wiederholung des Verfahrens wird man die Gleichung fort- 
gesetzt herabdrûcken, bis schliefilich aile Eoeffizienten bestimmt 
sind, was zu leisten war. 

Fûnfte Aufstellung. Lehrsatz. 

Zwischen beliebig vielen gegebenen GrôBen a, h, 
c, d, . . ., unter denen es keine gleichen gibt, und ebenso 
vielen gegebenen GrôBen irgendwelcher Art a, (i,y,S, . . . 
mag eine beliebige Beziehung, von jeder zu jeder, be- 
stehen. Nennt man unter den ersteren GrôBen irgend- 
eine unbestimmt x und die ihr entsprechende unter 
den letzteren unbestimmt y, so wird dièse Beziehung 
sich immer durch eine Gleichung von folgender Form 
ausdriicken lassen 

y =z f-\^ gx + hx^ + kx^ + ... 
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Sie bat auf der anderen, dem y gleichen Seite so viele 
Glieder, als gegebene GrôBen a, b, c, d, . . . oder a, /S 
y, Sj . . . àa, sind. 

Es seien drei Grôfien a, b, e gegeben und ebenso drei an- 
dere u, ^, y. Icb bebaupte, dafi sicb Eoefûzientenwerte f, g, h 
finden lasseD, die der aufgestellten Bebauptung gentigen. Mau 
bat nâmlicb zur Âuffindung dieser drei Werte drei verscbiedene 
Gleicbnngen 

cc = f + ga + ha^, 

p^f + gb + hb\ 

y = f+ 90 +hc^. 

Daber lassen sicb f, g, h nacb bekannten Gesetzen der 
Ânalysis finden. Nacb âbnlicber Metbode wird sicb die Rich- 
tigkeit der Aufstellung ergeben^ wenn irgendeine andere Anzahl 
von GrôBen .vorliegt, was zu beweisen war. 

Sechste Anfstellang. Problem. 

Die krumme Linie paraboliscber Natur zu finden, 
die durcb beliebig viele gegebene Punkte bindurcbgebt. 

Es seien jene gegebenen Punkte -4, J5, C, D, . . . , und man 
ziebe von ibnen nacb einer beliebigen geraden Linie EL die 
Geraden AHj BJ, CK, DL in einer solcben beliebig gegebenen 




Fig. 4. 

Bicbtung, dafi nur eine einzige von diesen denselben Punkt 
jener triflft Man denke sicb AH, BJ, CK, DL als Ordinaten, 
die auf den Abszissen OH, OJ, OK, OL steben, welcbe von 
einem gegebenen Punkte O aus gemessen werden. Irgendeine 
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von jenen werde unbestimmt y genannt, irgendeine von diesen x. 
Es ist dann nach der vorigen Aufstellung klar, dafi die Be- 
ziehung zwischen Ordinale und Abszisse, von welcher Art aie 
auch sein mag, immer durch folgende allgemeine Gleichong 
ansgedrûckt werden kann 

y =s f^ gx + hx* + kx^ + ... 

Hier bezeichnen f, g, h, k, ... unverânderliche GrôBen, die aber 
noch nicht bestimmt sind. Da nun jene Gleichung sich auf eine 
Kurve parabolischer Natur bezieht, so steht fest, daB eine solche 
Kurve immer durch gegebene Punkte gelegt werden kann. 

Erste Lôsnng. 

Wenn die Lange irgendeiner neuen Ordinale PM gesucht 
wirdy die auf der gegebenen Abszisse OP sleht, so schreibe 
man a, /?, y, S beziiglich flir — HA, + JB, — KC, + LD und 
ebenso a, ô, c, d fur — QH, + OJ, + QK, + OL und z flir 
+ OP. Dann wird das nach der 3. Aufslellung zu ermitlelnde 
^ gleich der gesuchlen Ordinale PM sein, was zu leislen war. 

Zweite Lôsnng. 

Wenn die Gleichung zu der Kurve gewtinscht wird, so lâBt 
sich auch unter denselben Vorausselzungen nach der 4. Aufstel- 
lung leicht beslimmen, was zu leislen war. 

Zusatz. 

Wenn demnach beliebige GrôBen und ebenso viele andere, 
von denen sie abhângen, gegeben sind, so wird sich irgend- 
eine gewûnschte passend inlerpolieren lassen. Man braucht 
nur jene gegebenen GrôBen durch Ordinaten AH, BJ, CK, 
DL, ... darzustellen, die anderen aber, von denen die iibrigen 
abhângen, durch Abszissen OH, OJ, OK, OL, . . . und die ein- 
zuschiebende GrôBe durch die Ordinale PM, 

Es wird uns nicht gereuen, ein BeispieP fur dièse Mé- 
thode beizuftigen. Von der GrôBe 1 + seien die Potenzen 

1, 1 + 0, 1 + 20+0», 1 + 3 + 3 0' + OS . . . 
gegeben, sowie auch ihre Indizes 0, 1, 2, 3, . . . Zu finden sei 
von derselben GrôBe die Potenz 



(1 + 0)" 
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mit dem Index — . Man setze a, h, g, d, , . . beziiglich gleich 

den gegebenen Indizes 0, 1, 2, 3, ... und a, ^^y. S, ... gleich 
den gegebenen Potenzen 

1, 1 + 0, 1+20+ OS 1 + 30 + 30^+ OS ... 
Fur z schreibe man den Index — der gesuchten Potenz. Dann 

wird Ç (welches nach der 3. Aufstellung zu ermitteln ist) eben 
die Potenz sein, die gefunden werden solL Man erhâlt in der 
Tat, wenn man die Rechnung durchfiihrt, 
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^ - 1.2.3 ' 



Hieraus ergibt sich 

m 

= 1-1 Q-\ — — . Q2 H + . . . 

Siebente Aufstellang. Problem. 

Die Flâchen aller Kurven mit grôfitmôglicher An- 
nâhernng zn bestimmen. 

Man denke sich durch einige Punkte der vorgelegten Kurve 
die nach dem letzten Problem zu bestimmende parabolische 
Kurve gezogen. Ihre Flâche (die sich immer nach wohlbekann- 
ten Methoden erhalten lâfit), wird dann ziemlich genau mit 
derjenigen der vorgelegten Kurve iibereinstimmen. 

Es handle sich um die Kurve ABGD, Man beziehe auf 
die Gerade KR irgend wie viele Ordinaten AK, BL, CN, 
DR, . . . Aus diesen und den zugehôrigen Abszissen, vom 
Punkte K aus gemessen, ergebe sich die Gleichung 

y == f+ gx + hx^ + kx^ + ..., 

die die Natur der Kurve ziemlich genau ausdrûckt. Dann wird 
die der Abszisse KP oder x und der Ordinate PM oder y an- 
liegende Flâche AKPM 

^fx + \9x^ + \hx^ + \kx^+ ... 
sein. 
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Wiiuscht man lieber eine geometrische als élue arithme- 
tische Quadratur, so wird die folgende élégant genug sein.^ 

Man ziehe die Parallèle Ad zu der Geraden KR, Sie 
schneide die (nôtigenfalls verlângertén) Ordinaten LB, NG, RD 
in hj c, d. Man nehme LF, NO, RH so zu 5ô bzw. Ce, Dd, 
wie sich eine beliebige Lange KS [S nach Belieben gewàhlt) zu 
KL, KN, KR verhâlt, und trage LF, NG, RH von den Punk- 
ten L bzw. N, R in demselben Sinne ab, wie die Linien Bb, 
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Ce, Dd von den Punkten b, e, d ans liegen. Es gehe dann 
durch die Punkte F, O, H eine Linie und schneide ^^ in E. 
Man ziehe wieder Eh parallel zur Geraden KR und bilde unter 
Wiedérholung des Verfahrens eine andere Linie, die ^^ in / 
schneidet, dann eine weitere, die AK in O schneidet, bis man 
zuletzt auf eine Linie kommt, die gerade und parallel zu KR 
ist. Nun mâche man, daB KP, KS, KT, KQ in stetiger Pro- 
portion stehen und Ks die Hâlfte der Abszisse KP, Kt ein 
Drittel, Kq ein Viertel derselben ist. Man nehme die Ge- 
rade ka gleich der Geraden KA, ebenso ke zu KE, wie Ks zu 
KS, und ki zu KJ, wie Kt zu KT, endlich ko zu KO, wie Kq 
sich zu KQ verhâlt, usf. Dabei trage man ka, ke, ki, ko ein- 
zeln in demselben Sinne vom Punkte k aus ab, wie KA bzw. 
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KE, KJj KO vom Punkte K aus liegen. Nimmt man jetzt KX 
gleich der Samme aller nach derselben Seite vom Punkte k ans 
gerichteten Stticke oder der Differenz der entgegengesetzt ge- 
richteten, so wird das Rechteck KX* KP ziemlich genau gleich 
der Flâche AKPM sein, was zu leisten war. 

Erlaubt man^ dafi die Punkte 8 und P zusammenfallen, so 
wird die Eonstruktion etwas einfacher. Es wird dann nâmlich 
die Gerade k a wie vorher gleich der Geraden KA sein, die Ge- 
rade ke aber die Hâlfte der Geraden KE, die Gerade ki ein 
Drittel der Geraden KJ, die Gerade ko ein Viertel der Ge- 
raden KO usf. 

Anmerkung. 

Es bleiben noch gewisse andere Problème ûbrig, die einen 
nicht zu yeracbtenden Nutzen in der Astronomie haben. Von 
solcher Art ist die Bestimmung der grôfiten oder kleinsten Or- 
dinate, ebenso die Bestimmung der Schnittpunkte einer Kurve 
mit irgendeiner gegebenen Geraden. Dièses aber und anderes 
dieser Art ûbergehe ich hier, um ein Ubermafi zu vermeiden^ 
zumal es aus der bereits erschlossenen Quelle, deren Gebrauch 
sich liber die ganze Mathematik und Physik weithin ausbreitet, 
leicht abgeleitet werden kann. 

Nachschrift. 

Die obigen Aufstellungen, die im Jahre 1707 niedergeschrie- 
ben wurden, habe ich meinen akademischen Hôrern in einer 
ôflFentlichen Vorlesung im Jahre 1709 auseinandergesetzt Es 
war mir nâmlich unbekannt, dafi Newton eine Abhandlung ûber 
denselben Gegenstand verfafit batte ^, bis ich ein Ekemplar 
davon, gedruckt im Jahre 1711, von dem gelehrten Heraus- 
geber Herrn Jones als Geschenk erhielt, dem ich auch eine 
Abschrift meiner Entwicklungen durch seine Freundlichkeit be- 
wogen ûbersandte. Nach Durchsicht der Aufstellungen des be- 
rtihmten Mannes bemerkte ich, dafi ich meistenteils einen ver- 
schiedenen Weg gegangen war, was durch Vergleich seiner 
zweiten Aufetellung mit meiner vierten klar hervortreten wird. 
Besonders schôn und niitzlich erschien mir aber die Regel, 
nach der er aus vier gegebenen Ordinaten, die in gleichen In- 
tervallen liegen, die Flâche zwischen der ersten und letzten mifit. 
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Die ein wenig fortgesetzte Reihe solcher Regeln f&ge ich hier 
unten bei. Die erste Ton ihnen liefert die Flâche zwisdien den 
âoBersten von drei âquidistanten Ordinaten, die zweite gibt die 
Flâche zwischen den auBersten von vier âquidistanten und wird 
Newton yerdankt, nach der dritten findet man die Fl&che ans 
fiinf ëiinlich gegef)enen usw. Hierbei wird immer A ftir die 
Summe der ânfiersten Ordinaten gesetzt, d. h. fur die Somme 
der ersten und letzten, B fur die Summe der nâchstâuBersten, 
d. h. fiir die Summe der zweiten und vorletzten, G fur die Summe 
der folgenden, d. h. fiir die Summe der dritten und vorvorletzten. 
So wird die Reihe der Buchstaben D, E, Fy . . . fortgesetzt, bis 
man entweder zur mittelsten Ordinate gelangt, falls die Anzahl 
derselben ungeradë ist, oder zur Summe der beiden mittelsten, 
falls die Anzahl gerade ist. Endlich wird das Intervall zwischen 
den âufiersten Ordinaten, d. h. die Basis der zu quadrierenden 
Flâche, mit dem Buchstaben R bezeichnet®. 



Ordinaten- 
zahl 



Flâchen 



in 

IV 

V 

VI 

VII 

VIII 

IX 

X 

XI 



^V'-" 




'V'-" 




7^4 + 325+12(7 „ 




90 "^ 




19^ + 755+500 o 




288 '^ 




41^ + 2165+ 27 C + 2722) 
840 


'R 


751^1 + 35775+ 1323 0+29891) 



R 
17280 

989^+ 58885— 9280+ 10496Z)-4540J5; „ 

• ji 

l'8 350 
2857^+15 74 15+ 10800+ 193442) + 5778^ ^ 

89 600 
16067^ + 1 06 300 5 -48 525 0+272 4002) -260 5505 +4278682^ 

598 752 

VLSW, 



R 



Digitized by VjOOQIC 



26 Oatùss: Neue Méthode zur naherungsweisen Intégration 



Neue Méthode 
zur naherungsweisen Auffindung von Integralwerten. 

(Von Cari Friedrich Gauss, 1814.) 



1. 

Tinter den zur naherungsweisen numerischen Bestimmung 
Yon Integrâlen vorgeschlagenen Methoden behaupten einen aus- 
gezeichneten Platz die Regeln, welche, nach dem VorgaDge des 
grofien Newton, Cotes entwickelt bat. In der Tat, wenn man 
den Wert des Intégrais 

lydx, 



fi 



Yon x = g bis x ^ h genommen^ sucbt, so sind die Werte von y 
fur dièse âufiersten Werte des x und ftir beliebige andere in 
gleicben Zuwûcbsen fortscbreitende Zwiscben werte mit bestimm- 
ten numeriscben Koeffizienten zu multiplizieren. Dann wird das 
Aggregat der Produkte mal h — g das gesuchte Intégral mit 
um so grôBerer Genauigkeit liefem, je mehr Werte in dieser 
Opération benutzt werden. Da die Prinzipien dieser Méthode, 
welcbe die Geometer seltener als recht ist anzuwenden schei- 
nen, meiues Wissens nirgends ausfuhrlicher erlautert worden 
sind, 80 wird es nicht unangebracht sein, darûber etwas Yor- 
auszuschicken. 

2. 

Es sei n + 1 die Menge der Glieder, welche man zu be- 
nutzen beliebt bat, und wir wollen h — g =^ J setzen, so daB 
die Werte von x folgen.de sind 

9. 9+-, 9 + —-, 9 + ^^ ^8W. 

bis g + A, und denselben folgende Werte von y entsprechen 
A, A\ A", A'" usw. bis ^("h endlich werde unbestimmt 

X ^ g + At 

gesetzt, so daB y auch als Funktion von t betrachtet werden 
kann. Bezeichnen wir mit Y die folgende Funktion ^ 
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(nt- l)(nt- 2)(nt- 3) . . . (nt - n) 



+ A' 
+ A" 
+ A'" 



(-1) 

nt(nt 


1-2) 
-2)(nt 


(- 
-3) 


3) 


. (nt 


(-n) 
-n) 


1 (- 
nt(nt 


■l) (- 
-l)(nt- 


2) 
-S) 




■ (1 
.(nt- 


-n) 

-n) 


2 
nt(nt 


1 (- 
-l)(jit- 


1) 
-2) 




. (2 
.(nt- 


-n) 
-n) 


3 
nt(nt- 


2 1 

USW. 

-l)(nt- 


2). 




• (3- 
(nt- 


-n) 
« + 1) 



oder 



n(n - 1) (n - 2) . 

^^ j^(^) ^(a*) 



wo, wenn jit die einzelnen ganzen Zahlen 0, 1, 2, 3 ... n re- 
pràsentiert, 

TitA — ^ ^(^ ^ ~ 1) (^ ^ - 2) (^ ^ - ^) - - ' (y^ ^ ~ ^) 

~~ nt — Il * 

if 0") der Wert von T îûr nt = fi, 

Dann wird handgreiflich sein, daB Y eine ganze algebraische 
Funktion nter Ordnung von t bildet und ihre Werte fiir die ein- 

12 3 

zelnen n + 1 Werte von t, nâmlich 0, — ,—,— ... 1, den 

' ^ n^ n n 

Werten von y gleich sind. Weiter ist offenbar, daB, wenn 7' 
eine andere ganze Funktion bedeutet, die fiir dieselben Werte 
mit y iibereinstimmt, Y' — F fiir dieselben verschwindet und 

daher durch die Faktoren t. t ,/ , t ... < — 1 

teilbar ist und also auch durch deren Produkt (das von der 
Ordnung n + 1 ist). Daher muB zweifellos F', wenn es nicht 
geradezu mit Y identisch ist, eine hôhere Ordnung erreichen 
oder Y von allen ganzen die Ordnung n nicht ûberschreitenden 
Funktionen die einzige sein, die fiir jene w + 1 Werte mit y 
iibereinstimmt. Wenn also y, in eine nach Potenzen von / fort- 
schreitende Reihe entwickelt, vor dem Gliede, das <»» + * enthâlt^ 
iiberhaupt abbricht, so wird es mit Y identisch sein: wenn es 
aber wenigstens so rasch konvergiert, daB man die folgenden 
Glieder vemachlâssigen darf, so wird die Funktion Y inner- 
halb der Grenzen ^ = 0, < = 1 oder x ^ g, x = h die Stelle 
von y vertreten kônnen. 
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8. 

Nunmehr geht unser Intégral 

[ydx liber in A jydt, 

genommen von ^ = bis t = 1, wofiir wir nach dem soeben 
Bemerkten 

jfrdi 
einsetzen werden. Durch Entwicklung also von 7^) in 

wird 

jT(f*)di, 
von t = bis / = 1 , 

== h 1- USW. 

sein. Wenn man dièse GrôBe = M^'^ B^^ setzt, so wird das 
gesuchte Intégral 

= A{AR+A'R'+ Â'R" + A"'R'" + usw. + ^(«) ^(~)) 
sein. 

Wir woUen beispielsweise die Berechnung des Koeffizien- 
ten R'' flir n = 5 hinzufiigen. Es ergibt sich hier 

r'= 55^5- 13.5*/*+ 59.5»<3_ i07.5»/2+ 60-5./, 

if"=2x 1 X(-l)X(-2)x(-3) = -12. 

Daher 

-12iî"=iHi_i625 + i5!i_^5!^ + 150=-ii, 

also 

p// 25 

Die Berechnung wird etwas kûrzer, wenn man 2 < — 1 = w 
setzt Dann ergibt sich 

~~ 2"(nw + n--2/w) 

Setzen wir 

(n*^' ~ n«) » [n't^« - (n - 2)*] . [n^u* - (w -- 4)»] > [n^u^ - (n - 6)«] . . . _ ^, . 
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wo der Zâhler aufhôren muB mit . . . (w^w* — 9)(n*tt' — 1), weim 
n ungerade ist, oder mit . . . (n'w* — 4)nw, wenn n gerade ist, 
und es wird 

2" 

Nonmehr ist das Intégral 

von / = bis ^ s 1 genommen, gleich dem Intégral 

von tt » — 1 bis tt = + 1. 
Wird daher 

U^) = aw«-l + ^W«-3 + yt^n-6 + du^-1 + U8W. 

gesetzt (von selbst ist nâmlich klar, dafi die Potenzen u^"^, 
w*-*, w«-6 usw. fehlen), so wird von dem Intégral der Teil 

J 2'* + l 
fiir einen ungeraden Wert von n verschwinden, der andere Teil 

r i2fi--n)U^U u 
J 2" + ^ 

dagegen fur einen geraden Wert. Also wird das Intégral 

fT(f*)dt 
fur gerades n , 

ftlr ungeradea » aber 

In anserem Beispiel bat man 
U" = (25 «» - 25) (25 tt* - 9) =. 625 m* - 850 «» + 225 
nnd demnacb 

_12iî" = -±-(l25-i|«H-225)— f-, 

wie oben. Es ist ntitzlich zu beachten, dafi 
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wird, und demnach 

jT(^'f*)dt^±fT(f*)dt, 

wobei das obère Zeichen fur gerades n, das untere fîir uDgerades 
gilt. Da man leicht einsieht, daB ebenso 

so wird stets 

sein oder von den Koeffizienten 

JB, R\ B", R'" . . . i2(«) 

der letzte dem ersten gleich, der vorletzte dem zweiten und so 
weiter fort. 

4. 

Die von Cotes bis zu n = 10 berechneten numerischen 
Werte dieser Koeffizienten schreiben wir . aus der Harmonia 
Mensurarum hierher. 

Fiir 7i = 1 oder zwei Glieder 

Fur w = 2 oder drei Glieder 
R-K -~, R =- . 
Fur w = 3 oder vier Glieder 

Fiir w = 4 oder'fûnf Glieder 

^=^ -W ^=^ =-45-' ^=T5-' 
Fur w = 5 oder sechs Glieder 

R^^R'f^.}^^ R' = R''''=—, ^"«i?'"=-?^. 

288 ' 96 ' '144 

Fiir w = 6 oder sieben Glieder 

840 ' 35 ' 280 ' 
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Fur n = 7 oder acht Glieder 

^""^ "17 280' ^-^ ""17280' ^ -^ ~ "âïÔ"' 
p/// P"" ^ 2989 



17 280 

Fiir n = 8 oder neun Glieder 



28350 ' 14175 ' -c* — ^^^^^ , 

Ji'-^JIY^ 5248 ^iv^_ 454 



14175 ' 2835 

Fût w = 9 oder zehn Glieder 



«9 600 ' 89 600 ' 2240 ' 

^'" = jRVI =^ i???_ ^ 7JIV ^ jijjv , 2889 



5600 ' ^ 44800 

Fiir w = 10 oder elf Glieder 

. ""^ - 5Qg^52' ^ -^ -149688' ^ ""^ 199581' 



1:^474 ' 11088 ' 24948 

5. 

Da die Formel 

A [AB + a:R' + A'R" + A'"R"' + usw. + A^^)R^^)) 
das Intégral 

Jydx 

von x^ g bis a; = ^ + J oder das Intégral 

àÇydt 

von < = bis < = 1 zwar genau liefert, so oft y, in eine Reihe 
entwickelt, die Potenz <*» nicht iiberschreitet, aber nur nâhe- 
rungsweise, so oft y dariiber hinausgeht, bleibt tibrig, daB wir 
den Fehler, den die nâchstfolgenden Glieder hineinbringen, be- 
stiinmen lehren. Bezeichnen wir allgemein mit Ai^"») die Diffe- 
renz zwischen dem wahren Wert des Intégrais 



Çt^^dt, 
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von < = bis < = 1, und dem Wert, der aus der Formel her- 
Yorgelit, 80 dafi 

ifc = 1 -i2-i?'-.^"_i?'"-U8W. -i?('»), 

k' =-l--l(72'+2/?"+ 3^'"+u8w. +7Ï i?(«)), 

r = -1 - -L [E' + 4i?" + 9^'" + U8W.- + n»i?(«)), 

r' = i- - ^(i2' + 8JS" + 21 K" + U8W. »»^('»))* 

U8W. Es ist daher offenbar, daB, wenn man y in die Reihe 

K + K't + K"t^+K'"t^+ U8W. 
entwickelt, die Dififerenz zwischen dem Wert des Intégrais 

\dt 

und dem angenâherten Wert der Formel ausgedrtickt wird darch 

Kk + K'K + K''k" + K'''r' + usw. 

Aber handgreiflicherweise werden A;, k\ k" usw. bis A;^*») von 
selbst = 0: Die Korrektion der angenâherten Formel wird 
also sein 

Die Beschafifenheit der GrôBen /:('*+!), A;('»+2) usw. werden 
wir unten genauer erforschen; hier môge es gentigen^ die nu- 
merischen Werte der ersten oder zweiten, fur einzelne Werte 
von n, beigefiigt zu haben, damit danach der Grad der Ge- 
nauigkeity den die angenâherte Formel bietet, geschâtzt wer- 
den kann. 

Fiir n = 1 haben wir 

Fur » = 2 finden wir 

120"' '*' ^""lê" 



1/" — fr"" ■— — ^ frv __ _ ^ 



Fur n « 8 wird 



t"" — — ^ tv « __ ^ 
'^ "■ 270' '^ " 108 



* In der Ausgabe der WW. Gôttingen, 1866, steht fîllsclilioh 
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Fur n = 4 wird 
Fiir n = 5 wird 



'^ - "' '^ - 2688 ' '^ " 768 



/g VI -^ . . ^^ » jfcvn -- _ . ^^ 



52500 ' 15000 

Fur n = 6 wird 

'^ -^, '^ - 38880 . ^ - 8640 

Fût w = 7 wird 

JkVIII -_ _ ^^"^ JfclX ^ _ ^^"^ . 

10588410' 2352980 

Fiir w = 8 wird 

A;IX==0, ^^=- 17301504**' ^^^ "" "" 3145728 

Fiir n = 9 wird 

1.x -- _ 865 I.XI ^ _ 865 

"' RQ1 Qf^l onft > 



631351908' 114791256 

Fiir w = 10 wird 

JUXI-.0 7,^TT _ 26927 JkXIII^_ 26927 



136500000000 ' "" 21000000000 

Fiir einen geraden Wert von n sehen wir hier iiberall 

A;('» + l) = 

werden und auBerdem 

Jk(n+3)^'i?_±_iA;(n+2). 

fiir einen ungeraden Wert von n aber ergibt sich iiberall 

jfc(n+2) ^ ^L±_ijfc(n+1). 
2 

Der Grund dieser Erscheinungen ist leicht ans den folgenden 
Ûberlegungen zu entnehmen. 

Bezeichnen wir allgemein mit /("*> die DiflFerenz zwischen 
dem wahren Wert des Intégrais 



Sk-^)'"' 



* In der Ausgabe der WW. Glôttingen, 1866, steht fôlschlich 
In Nr. 6 hat aie wieder den richtigen Wert A;^^. 



52300 

37 
** In der genannten Ausgabe steht wiederum falschlicb — ■ ■ . 

Kowale wBki, Newton luw. 3 
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yon t = bis ^ = 1, und dem Wert, den die angenàherte For- 
mel darbietet^ so daB man bat 

' -;('-i^-[(-i^+(i-ir«'+(i-Tr«" 

das Intégral yon < = bis ^ = 1 genommen. 

Oflfenbar wird fur einen ungeraden Wert von m sowohl 
der wahre Wert des Intégrais als auch der angenàherte Wert 
verschwinden: es wird also 

r=0, Z"'=0, /v = 0, /vii^o usw. 

oder allgemein /("») = flir einen ungeraden Wert von m. Fur 
einen geraden Wert von m aber geben wir der Formel solche 
Gestalt 

+ usw. +2«i?(^''"')+i?(^''"0], 
wenn gleichzeitig n gerade ist; oder dièse 
lim) = _L L-L-^ ^ln^B + (n^ 2>rR' + in - 4)«»i?" + usw. 

wenn gleichzeitig n ungerade ist. 

Wenn daher durch Entwicklung von y in eine nach Po- 
tenzen nach t — \ fortschreitende Reihe herauskommt 

2/ = L + r(i - X) + L"{t - i)2+ r^{t - i)3+ usw., 

80 wird bei dem angenâherten Wert des Intégrais iyàt, von 
^ = bis / = 1, die Korrektion anzubringen sein 
Ll + L"r+ L"" V" + LVi /vi + usw. 

oder vielmehr wird, da ^'") fiir jeden ganzzahligen Wert, der 
nicht grôBer als n ist, notwendig verschwindet, die Korrektion 
im Falle eines geraden n 

J^(n+2)/(n+2) ^ 2/('»+*) /(»+*> + L(" + ^^)^"+Vï) + USW. 

laaten oder 
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im Falle eines ungeraden n. 

Sehr leicht lassen sich nun die Korrektionen Z^«) auf A;^*») 
zuruckfiihren und umgekehrt. Da nâmlich 

80 wird 

A 4 1 • ^ 

sein. 

Und ebenso ergibt sich 

Aus letzterer Formel werden die Glieder herausfallen, bei denen 
l mit ungeradem Index behaftet ist: jede von beiden ist aber 
nur bis zum Index n+1 (einschlieBlich) fortzusetzen. Oflfenbar 
werden wir also haben: 
Pur gerades n 

A;(n+1) = 0, 
2 

Fur ungerades n 

A;("+l) = ^«+1), 
2 

woraus die oben angezeigten Beobachtungen entspringen. 

' 6. 

Im allgemeinen ist es darum vorteilhafter, bei der Anwen- 
dung der Cotesianischen Méthode dem n einen geraden Wert 
zu geben oder eine ungerade Menge von Gliedem zu benutzen. 
Gar wenig wird nàmlich die Genauigkeit yermehrt werden, wenn 
wir statt des geraden Wertes von n zu dem nâchst grôBeren 
ungeraden aufsteigen, da der Fehler von derselben Ordnung 
bleibt^ mag er auch mit einem etwas kleineren Koeffizienten be- 
haftet sein. Dagegen wird durch Fortschritt von einem un- 

3* 
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geraden Wert des 7i zum nâchstfolgenden geraden der Fehler 
um zwei Ordnungsstufen fortgeschoben werden und ûberdies 
wird der merklicher verminderte Koeffizient die Genauigkeit 
erhôhen. So lâBt sich, wenn fiinf Glieder zur Anwendung 
kommen, oder fiir n = 4, der Fehler nàherungsweise durch 

7^^^^^ ausdriicken oder durch — i^^^^i^^®^; wenn man 

w = 5 setzt, 80 wird der Fehler nàherungsweise — K^^^ 

^^®^ "■ t;op;nn ^^^^ ^®^^ ^^^ ^^ nicht einmal bis zur Hâlfte des 
friiheren herabgedrûckt; bei n = 6 dagegen wird der Fehler 

annâhemd =~ ggiâô"^^^' ^^®^ ~ sskô"^^^^' ^^ ^^® Genauig- 
keit wird um so grôBer, je rascher die Reihe, in welche man 
die Funktion entwickelt, schon von selbst konvergiert. 

7. 

Nachdem wir dies liber die Méthode von Cotes voraus- 
geschickt haben, schreiten wir zur allgemeinen Erôrterung, wo- 
bei die Bedingung, da6 die Werte von x eine arithmetische 
Reihe bilden, fortfâllt. Wir nehmen also das Problem in An- 

grifif, den Wert des Intégrais jy dx zwischen gegebenen Grenzen 

ans irgendwelchen gegebenen Werten von y entweder genau oder 
nàherungsweise zu bestimmen. Fùgen wir hinzu, daB das In- 
tégral von X =^ g bis a; = ^ + J zu nehmen ist, und fuhren wir 
statt X eine andere Variable 

ein, so daB es notwendig ist, das Intégral 

A^ydt 

von / = bis ^ = 1 zu untersuchen. Es môgen n + 1 ge- 
gebene y -Werte 

A, A', A\ A"' . . . ^^"' 
den ungleichen ^-Werten 



entsprechen, und wir wollen mit F folgende ganze algebraische 
J'iinktion nter Ordnung bezeichnen*: 
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^ (t-a')(t-a"){t-a"') .. 


. (^ - a(»>) 


(a — a') (a — a") {a — a"') . . 
(t-a)(t-a''){t-an .. 


, . (a - a^">) 
. (t - a<»>) 


' "^ (a' - a) la' - a") (a' ~ a'") . . 




' ^ (a" -a) (a" -a') (a" -a'") .. 
+ usw. 
. .(«) (/-«)(^-aO(^-a'0. 





Offenbar fallen die Werte dieser Funktion, wenn t irgend- 
einer von den GrôBen a, a, a" ... a^»»^ * gleichgesetzt wird, mit 
den entsprechenden Werten der Funktion y zusammen, woraus 
wir ebenso wie in Artikel 2 schliefien, daB Y mit y identisch 
ist, so oft y gleichfalls eine ganze algebraische Funktion ist, die 
die Ordnung n nicht iibersteigt, oder wenigstens ftir y eintreten 
kann, wenn y in eine nach Potenzen von t fortschreitende Reihe 
verwaudelt eine so starke Eonvergenz darbietet, dafi man die 
Glieder der hôheren Ordnungen vernachlâssigen darf. 

8. 

Nunmehr woUen wir zur Ermittlung des Intégrais I Ydt 
die einzelnen Teile von Y betrachten. Bezeichnen wir das 
Produkt 

{t - a){t - a)(t - a''){t - a'") ...(<- a^«>) 

mit T, und durch Entwicklung dièses Produkts werde 

T= ^«+1 + £zr + a'<'»"^ + a"^'»"^ ^ usw. + a^«^ 

Der Zàhler des Bruches, mit dem A in dem ersten Teile 

T 
von Y multipliziert ist, wird = ; die Zâhler in den fol- 

genden Teilen sind ebenso 

T_ T T 

t-a'" t'-a"' t^a'" ^®^" 

Die Nenner aber sind nichts anderes als bestimmte Werte 
dieser Zâhler, wenn beziehungsweise 

^ = a, / — a', < = a" y t = a'" usw. 



* Hier hat die Gôttinger Ausgabe von 1866 wieder einen Druck- 
fehler (a^w) statt a^«)j. 
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gesetzt wird: wir wollen dièse Nenner beziehungsweise mit 3/, 
M'j M'\ M'" usw. bezeichnen, so da6 
^ AT , AT , A"T , , -4^7 



Mif - a) ^ if' (^ - a') ' M" (t - a") •*•"•• Jtf («) (^ - a<»0 
ist. Da r = wird fur t = a, haben wir die identische Gleichung 
a"+i + cca^ + 6f'a«"^ + of"a«-"2 + usw. + f^«> = 

and demnach 

Hieraus wird durch Division mit t — a 

—^ = <« + a r~^ + a^ r-2 + a» /»-» -f- usw. + a« 

+ a f'"'^ + aa t""-^ + aa^ ^""'^ + usw. + u a""^ 
+ a' r-2 4- a' a ^«-3 + usw. + a' «""^ 
+ a" ^""3 + usw. + a" a*»-^ 
+ usw. usw. 

+ a^"-^^ . 
Der Wert dieser Funktion fur / = a betrâgt 
= (n+l)a" + wûfa""^ + (n — l)c^'a«--*+(w— 2)c^"a"^'+usw.+ c^^~"^^ 

j m 

Daher ist M gleich dem Werte von -rr- fur ^ = a, wie aueh 

et* 

anderweitig^ feststeht. Ebenso werden M\ M", if'", usw. Werte 

dT 
von -jr- sein fiir t ^ a!, t = a", ^ = a" usw. 

, von 

t — d 

< == bis ^ = 1 

+ ^ + ::^ — r + ::. -, + "sw. + o» 



« + 1 w w— 1 w — 2 



'nn— lw~2 

H + USW. + a a"" ^ 

' n — 1 n — 2 

+ -^^ + usw. + «"a«-3 

n — 2 



+ USW. USW. 
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Dièse Glieder werden wir foIgendermaBen umordoen: 

o» + f^a«-i + ûf'a«-2 + ^"a«-3 + usw. + r^^"-'> 

3 

+ -^(a"~' + a a"-* + a a"-' + usw. + «'""*') 
+ usw. 

+ ^ 



n+ 1 



Offenbar kommt aber dieselbe GrôBe heraus, wenn in dem 
Produkt aus der Multiplikation der Funktion T mit der unend- 
lichen Reihe 

nach Ausschaltung aller Glieder, die die Potenzen von t mit 
negativen Exponenten enthalten (oder kûrzer, in dem Teil des 
Produkts, der eine ganze Funktion yon < ist), / fur a geschrieben 
wird. Nehmen wir also an, es werde* 

T{t-' + Y^""' + y ^"' + T^"* + ^s^-) = ^'+ ^"' 

80 daB T die in dem Produkt enthaltene ganze Funktion von t 
ist, r" hingegen der andere Teil, nâmlich eine absteigende und 
ins Unendliche laufende Reihe. Danach wird der Wert des 

Intégrais f-f^ , von < = bis / = 1, dem Wert der Funktion T 

gleich sein fQr t = a. Wenn wir daher die bestimmten Werte 

der Funktion 

r 



m 



* „E8 wird kaum nôtig sein zu bemerken, dafi die Zeichen T, T', T" 
hier in andçrem Sinne genommen werden, als in Artikel 2." 
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fur ^ = a, ^ = a', < = a", t = a'" usw. bis t = a^»^ beziehuDgsweise 
mit R, R\ R'% i?'". . . . i2^*»^ bezeichnen, so wird das Intégral 
^Ydi, yon < = bis < = 1 , 

^RA + B'A'+ R!'Â'+ usw. + R^ J^^ , 
was mit A multipliziert den wahren oder den angenâherten 
Wert des \y dx, von x =^ g bis x ^ y + A, liefern wird. 

9. 

Dièse Operationen werden wesentlich leichter zustande ge- 
bracht, wenn man statt der unbestimmten GrôBe t eine andere u = 
2^—1 einfùhrt. Wir schreiben auch zur Abkûrzung 6 = 2a — 1, 
è'= 2a — 1, 6"== 2a"— 1 usw. T môge, nach Einsetzung des 

Wertes y ^ "^ T ^^^ ^* ^^^^8®^®^ i^ "ïhT ^^®^ ®^ ^®^ 

Z7= (w - 6)(w - b'){u - 6")....(w - b^). 

Es wird daher 

(^r ^ 1 dv 

dt 2" rftt 

, sein, und somit werden M, M\ M" usw. bestimmte Werte von 

_1_ dU 

wenn man der Reihe nach w = 6, u = b', u = b" usw. setzt. 
Da die Reihe 

nichts anderes ist als 

SO wird sie durch die Einsetzung t^—u-\- — notwendig ûber- 
gehen in 

2w-i + \u-^ ^ ^vT^ -\' ^u-^ usw. 

Wenn wir daher setzen 
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derart, daB U' die in diesem Produkt enthaltene ganze Funktion 
von u ist, U" hingegen der andere Teil, nâmlich eine ab- 
steigende nnendliche Beihe, so ist offenbar 

r' aber wird handgreiflicherweise als ganze Funktion von t 

durch die Einsetznng < = -— w + — - notwendig eine ganze Funk- 

tion von u hervorbringen: T' dagegen, welches nur négative 
Potenzen von t enthâlt, vnrd durch dieselbe Einsetznng lediglich 
négative Potenzen von u erzeugen. Deswegen wird U' nichts 
anderes sein als das durch dièse Einsetznng umgeformte 2'^T' 
und ebenso £7" ans 2'»T" gewonnen werden. Es wird also kein 
Unterschied sein, ob wir in 



t = a einsetzen oder in 






u ^ b machen, Daraus entnehmen wir, daB auch R, R\ R'\ 
R'" usw. bestimmte Werte der Funktion 

("Sir) 
fur w = 6, w = ô', w = V\ u = V" usw. sind. 



10. 

Bevor wir weiter fortschreiten, woUen wir dièse Regeln 
durch ein Beispiel erlâutem. Es sei n == 5, und setzen wir 

fest Daraus ergibt sich 

Durch Multiplikation mit 

r 1 + -1 r ^ + y r « + ~ <-* + usw. 
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erhalten wir 



r=t^-^t* + ^t^-^t^ + 4^t- '' 



l '80 20 • 7500 7Ô00 

Die Werte der Koeffizienten E, It\ R", R"\ K"\ R'"" werden 
demnach ausgedrûckt durch die gebrochene Funktion 



5 80 20 7500 7500 

^5 5 ^625 625 

wo ftir t der Reihe nach die Werte 0, —, —, -r-» -r» 1 ^i"- 



zusetzen sind. Viel kiirzer ist die andere Méthode, welche 
liefert 

5 ^ 625 625 

Ti' S 16 « , 277 

wonach i?, R\ R" usw. die Werte der gebrochenen Funktion 

4 16 . 277 



15 ÎH75 



^ ^ 28 , 518 
5 625 

3 1 

fûrw = — 1, w = — -, u = — - usw. Beide Methoden liefem 



dieselben Zahlen, welche wir in Artikel 4 aus der Harmonia 
mensurarum angef&hrt haben. Ubrigens werden in einem Falle, 
wie ihn dièses Beispiel darstellt, wo a, a', a" usw. rationale 

jtrp 

GrrôBen sind, die Werte des Nenners -jr- bequemer in der ur- 

spriinglichen Form berechnet, nàmlich 

(a - «^(a - a")[a - a'") . . .(a - é"^) 

fur t ^ a und ebenso im ûbrigen. Dasselbe gilt vom Nenner 

-7— , der fur w = 6 wird 
du ' 

« (6 - h'){h - b")[b - 0---(^ - ^^"^ 

11. 

So oft a, a', a" usw. entweder teilweise oder ganz irrational 
sind, wird die Umformung der gebrochenen Funktion, aus der 
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wir die Zahlen R, R\ R" usw. ableiten, in eine ganze Funktion 
nùtzlich sein. Da die Prinzipien solcher Umformung in den 
algebraischen Biichern nicht zu iinden sind^ wollen wir sie an 
dieser Stelle kurz aaseinandersetzen. Sind drei ganze Funk- 

tionen Z, f, Ç der nnbestimmten GrôBe x vorgelegt, so wird 

z 

also eine ganze Funktion gesucht, welche die gebrochene -^ 

vertreten kann, insofern fiir x irgendeine Wurzel der Gleichung 
^=0 angenommen wird. Wir wollen aber voraussetzen, da6 K 
fiir keinen dieser Werte von x verschwinde oder, was auf das- 
selbe hinauskommt, da6 f und Ç keinen gemeinsamen uii- 
bestimmten Teiler enthalten. Die Exponenten der hôchsten 
Potenzen von x m ^ und ^' werden wir mit k, h' bezeichnen. 
Es werde in gewohnter Weise f durch Ç dividiert, bis 
die Ordnung des Bestes unter k' berabgedrtickt ist; der Rest 
werde = -^ gesetzt und seine Ordnung = k'\ so daB y^^" 
das hôchste Griied des Restes ist; den Quotienten der Division 
wollen wir = y setzen. Ebenso môge aus der Division der 

Funktion ^ durch Ç' der Rest -^ von der Ordnung k"\ der 

Quotient ^ hervorgehen; sodann wiederum aus der Division 

der Funktion Ç' durch Ç" der Rest -^ von der Ordnung k"" 

und der Quotient ^ und so fort, bis man in der Reihe der 

Funktionen ^', f", Ç'" usw., welche einzeln ihr hôchstes Griied 
mit dem Koefiizienten 1 behaftet haben werden, zu Ç^ ^ \ 
gelangt. DaB dies endlich eintreten muB, lâBt sich leicht ein- 
sehen, da jede beliebige von den Funktionen f, ^, ^', ^" usw. 
mit der voraufgehenden keinen gemeinsamen nnbestimmten Teiler 
haben kann und daher gewiB keine Division ohne Rest aus- 
fiihrbar ist, Solange der Divisor eine hôhere Ordnung als hat. 
Wir werden also die folgende Reihe von Gleichungen erhalten 

r =A ? -;> r , 
r =^' r -/ r , 

usw. bis 

nm) «- ;[(to— 2) «m— 2) ^dn— 2)«m— 1) 
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wo C\ Ç'\ r" *"Ç'^ ganze Funktionen des x yon der Ordnung 
r, r; r"...ifc^"'>; die ZahleD A;', r, r'...A^«> bestândig ab- 
nehmend bis zur letzten U"^ ^Q\ p, p\ p\ /" usw. gleichfalls 
ganze Funktionen des x von der Ordnung k — k\ A;'— k'\ k"— k'", 
k'"— k"" usw. (mit Ausnahme des PaUes k < k', wo oflfenbar 
jo = gesetzt werden mu6). 

Nach solchen Vorbereitungen bilden wir eine zweite Eeihe 
yon ganzen Funktionen des x, nâmlich ti, ?/, rf\ rf" . . . rj^^'K Und 
zwar woUen wir rj = 1, i?'= festsetzen, die ûbrigen aber leiten 
wir einzeln aus den beiden voraufgehenden nach demselben 
Gesetz ab, durch das die Funktionen f, ^, ^\ ^" usw. unter- 
einander verkntipft sind, nâmlich nach den Gleichungen 

f]'' ^X f] --p ri' , 

fj := ?. 7/ —p 7j , 
lin -ïii ^11 ^11 III 



'//// •\in III III iiii 



USW. bis 

fjim) 5» ;^(m— 2) ^ (m— 2) _ «j(ni— 2) ^(m— 1) ^ 

Oflfenbar ist hier V'=^ ^on der Ordnung 0; f/"=--A;7' von 
dèr Ordnung Aï'— A;" und ebenso die folgenden 17"", ly'"" usw. 
beziehungsweise von der Ordnung A;'— A;'", A;'— A;"" usw., so daB 
das letzte t/*^ zur Ordnung A;'— A;^"»~^* aufsteigt 

Ferner wollen wir eine dritte Eeihe von Funktionen be- 
trachten, f-^17, T- ^^'» T-f'A T'-fV" iisw., wo unter 
drei beliebigen aufeinanderfolgenden Gliedem offenbar eine âhn- 
liche Relation bestehen wird, n'dmlich 

r -?v' =A ic -cv)-P (r -?^'), 
r - f ^'" = A' (r - ?^') -/ (r - ?^"), 
r'- ?^""= r (r- ?^") - ^"(r - so- 

Nun wird die erste dieser Funktionen =s 0, die zweite = Ç: 
daraus lâBt sich leicht schlieBen, daB die einzelnen durch Ç 
teilbar sein werden. 

Hieraus folgt ohne Scbwierigkeit, daB an Stelle des Bruches 
— die ganze Funktion Zt/^"»^ angenommen werden kann, insoweit 

wenigstens dem x nicht andere Werte erteilt werden als die, 
welche Wurzeln der Gleichung f» sind: oflTenbar muB nâm- 
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lich die Differenz — ^ — . '^ fiir solchen Wert von x ver- 
schwinden, da 1 — ^rj^""^ = ^"»^ — f i?^™^ durch ^ teilbar ist 

An Stelle der Funktion Z?/*"^ wird auch deren Rest aus 
der Division mit ^ treten kônnen, der von niedrigerer Ordnunng 
sein wird, als die Funktion ^. 

Ubrigens lâBt sich dieser Rest bequemer durch folgenden 
Âlgorithmus unmittelbar ausfindig machen. Es môgen fol- 
gende Gleichungen gebildet werden 



z 


= ?' r +z' , 


Z' 


= q" r + Z" . 


Z" 


= q'" r" + Z'" , 


Z'" 


^ q"" t"" -X- Z"" , 



usw. bis 

indem man nâmlich nach der Reihe Z durch f ' dividiert, darauf 
den Rest der ersten Division Z' durch Ç", dann den Rest der 
zweiten Division durch f '" und so fort. Da der Rest stets zu 
einer niederen Ordnung gehôrt, als der Divisor, wird die Ord- 
nung der Funktionen Z\ Z'\ Z"\ Z"" usw. beziehungsweise 
niedriger sein als A;', li\ k"\ k"" usw.; die letzte voUends Z^»»^ 
wird notwendig = 0, da der Divisor Ç^"*^ = 1 ist Wir haben also 

Insofern aber fur % nur Wurzeln der Gleichung Ç' = an- 
genommen werden, wird 

wonach unter derselben Beschrânkung 

— - = g^-f-ç ^ -\. q ff -j- USW. + q^^ fj^^ 

herauskommt Die Ordnung aber dièses Ausdrucks wird not- 
wendig unterhalb k' sein: da nâmlich die Ordnung der Quotienten 
g", q"', g"" usw. unterhalb A;'- r, r- F', r'~ F" usw. sein 
mu6, wird die Ordnung der einzelnen Teile q''fj'\ q'' r['\ 
<i"'ir usw. unterhalb A;'- A;", A;'- A;'", A;'- A;"" usw. sein. 

Endlich bemerken wir noch, wenn sich zufâllig unter den 

.Werten der Unbestimmten %^ die man in den Bruch — ein- 
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setzen muB, rationale mit iiTationalen gemischt finden, daB es 
zweckmâBiger sein wird^ jene von diesen zu sondern und dièse 
allein in der Gleichung ^'^0 zusammenznfassen. Fur rationale 
Werte nâmlich wird eine Recbnungserspamis nicht nôtig sein; 
fur irrationale aber wird die Rechnung desto einfacher sein, je 
geringer der Grad der ganzen Funktion ist, auf die man die 
gebrochene zurlickfiihren kann. 

13. 

Hier nun ein Beispiel fiir die im yorigen Artikel aus- 
einandergesetzte Umformung. Es sei die gebrochene Funktion 

50 ^ 283 , 256 
89 ^ 715 15015 



^ . 105 . 315 , 35 



Yorgelegt, in der z unbestimmt die Wurzeln der Gleichung 

, 21 ^5 105 3 35 ^ ^ 

* -IF* + m"* -729"*="" 
darstellt. 

Wenn wir hier aile sieben Wurzeln umfassen woUten, 

wûrden wir auf eine ganze Funktion sechster Ordnung kommen. 

Offenbar aber ist fiir den rationalen Wert ;?; = die Berechnung 

ne A 

des Bruches leicht und gibt den Betrag -r^^r-. Daher wollen 

wir unter Ausscheidung dieser Wurzel bei der Gleichung 
sechsten Grades verweilen: 

"^ 13 ^ + 143 "^ 429 " "• 

Dann sehen wir leicht voraus^ daB eine ganze Funktion vierter 
Ordnung herauskommen wird. Nunmehr ergibt sich aus der 
Anwendung der yorangegangenen Regeln folgendes: 

* 13 143 42» 

^ " " 13 * +T4F* 129"' 



?""= 1 . 



11 ' 83 

3 
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X 


13 
== l2-' 


r 


4719 
280 ' 


X" 


147 


8 ' 


V 


= 1. 


v' 


= 0, 


< 


13 

= 42 ' 



13 

, 4719 2 , 3333 

-^ " 280 ^ "^ 280 ' 

147 2 , '^'^'^ 

^ 8 ^ 88 



20449 « 14443 



^ 3920 ^ 3920 ' 

61347 . 127413 « . 120263 
' 640 1120 ^ 4480 

rr « 50 A . 283 o 256 , - 

^ = ^'--39r^-' + -7ïr^ -T5ôr^' ^ =^' 

' - ± 4 _ li 2 . 323 „ Jl_ 

■" 3 ^ 65 * "^ 5005 * ^ • ^ 3 ' 

2145 " "^ 45045 ' ^ 2145 ' 

** 3465 ' ^ "T 3465 ' 

Hieraus lâBt sich endlich die dem yorgelegten Bruch âquivalente 
ganze Fanktion ableiten: 

1859 . 1573 2 . '794'? 



16800 29400 . ' 39200 



18. 



Um den Grad der Genauigkeit zu bestimmen^ deren sich 
unsere Integralformel 

RA^ R'Jl + i2"^"+ usw. + ^<«>^<«> 

erfreut, woUen wir allgemein setzen 

Rù?^ + Rd"^ + i?V« + usw. + ^^»>a<'»^«» == — \— — A;^"»^ , 
» ' ' ' w + 1 ' 

80 daB ib^*^ die Differenz zwischen dem wabren and angenâherten 
Wert des Intégrais \i'^di ist, genommen von < = bis <= L 
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Wir werden daher nach Entwicklung der einzelnen BrQche ia 
Beihen erhalten 

= (1 - A) r-i + (1 - fcj r« + (i- - k") <-" 

= ri + i-r» + yr» + i r* + usw. - 0, 

wenn wir 

= A; r 1 + A;'r 2 + r r » + r'r* + usw. 

setzen oder vielmehr (da wir bereits wissen, dafi k, k\ k'\ k'" usw. 
bis U^ Yon selbst verschwinden mûssen) 

= ]k(n+l)^-(n+2) ^ ;j.(n+2)^-(n+3) ^ jj.(n+ 8) ^ - (n+ 4) ^ ^g^^ 

Durch Multiplikation mit T ergibt sich 

^ (t^^ + 7^ + -CT- + '^«-- + tS^) = ^+ ^"- ^®- 

Die linke Seite dieser Gleichung ist eine ganze Funktion des t 
von der Ordnung n, und ihre fiir < = a, t ^ a\ t ^ a" usw. be- 
stimmten Werte werden beziehungsweise MR, M'E% M"R'' usw.: 
deswegen mufi, da dasselbe von der Funktion T gilt^ wie aus 
der Bestimmungsweise der Zahlen R, R\ R" usw. ersichtlich 
ist, jene linke Seite der Gleichung notwendig mit' T identisch 
sein und somit r"= T0. Es entsteht daher aus der Ent- 

wicklung des Bruches -^ , und auf solche Art werden sich die 
Koeffizienten kf^'^^\ A;^»+2) ^g^ beliebig weit bestimmen lassen. 
Nach diesen Ermittelungen wird die Korrektion von unserem 

angenâherten Wert des Intégrais j ydt sein 

= k^n+l)X^n+l) ^ jfc(n+2) ^n+2) ^ ^g^^^ 

wenn die Reihe, in die y entwickelt wird, die Gestalt hat 
y==K+KU + K"fi + K"'i^ + usw. 
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14. 

Will man die Korrektion lieber durch Koeffizienten der 
nach Potenzen von t — -^ fortschreitenden Reihe 

y = L + L'{t--^)+L"{t- !)'+ L'" {t - -1)'+ usw. 

ausdriicken, so wird jene 

sein, wenn wir allgemein mit ï^"^ die Korrektion des angenàherten 
Wertes des Intégrais 



/( 



-i)'"' 



ausdriicken. Dièse Korrektionen l^"^ werden mit den Korrek- 
tionen Ic^"^ durch die Gleichung 

- T ÎT2T3 *^" ' + '^«^- 

verkntipft sein. Um jene aber unabhângig ermitteln zu kônnen, 
wollen wir erwâgen, daB die Funktion durch die Substitution 

ûbergeht in 

2 k (w~^ — w~^ + ti""' — M""* + usw.) 
+ ék' (m-2 — 2w-'+ 3w-"*- 4tt"-'* + usw.) 
+ 8 r (m-3 - 3 W-* + 6 M-» - 10 î*-« + usw.) 
+ 16 k'" (î*-* - 4 w-6 + 10 î^-« - 20 M-^ + usw.) 
+ usw. 
oder in 

2 AW-1+ 4 (ik'- i-Jî^^a + 8 (r- y • 2;k'+ T*)^"^ 

+ 16 (r'- i- . 8r+ i- . 3fc' - -i-t) W-* + usw. 
oder in 

oder endlich, da wir a priori wissen, daB /, /', T, T' usw. bis f^ 
von selbst verschwinden, in 

2«+2/(n + D|^-(n+2) ^ 2*» + ^ P"+® w'"^'»+* + 2*»+* ?'» + *^ tt"'^'*+** + USW. 
KowAlewskl, Newton usw. 4 
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Nun ist aber = -=r-; daher wird, weil T, T" durch die 

Substitution < = ^w + |- in ^^^ , — ^ iibergehen (Art. 9), die 

2 TJ" 
Funktion durch dieselbe Substitution in -yr- iibergehen. Wenn 

TJ" 
wir also die ans der Entwicklung des Bruches -jj- entstehende 

Reihe mit I> bezeichnen, so wird 

_[. 2«+3 /(«+3) ^ - (n+4) iig^^ 

sein, und auf dièse Art werden sich die Koeffizienten Z("+i>, 
Z('»+2) usw. beliebig weit ermitteln lassen. 

So finden wir in dem Beispiel des Art. 10 

rrn 116 . 304 _, 2576 -c 

(7 = w ^ w "* w * — usw. 

13125 28125 3U9 375 

^ 176 _, 832 _n 189 856 _,, 

fi = U ^ U ^ U ^^ — USW., 

13125 28125 4296875 ' 

und die Korrektion des angenâherten Wertes des Intégrais ist 
daher 

^ ^ LVi 1^ LViii ^?^? — Lx - usw. 



52 500 112500 137500000 

15. 

Der Koeffizient K^"^^ der in eine Reihe entwickelten Funk- 
tion y wird nach dem Satz von Taylor gleich dem Wert von 



1.2.3...W dC^ 1.2.3...7W dx"^ 

fiir / = oder a? = ^; ebenso ist der Koeffizient L("») der Wert 
desselben Ausdrucks fur t ^\ oder w = oder x = g + ^ J: 
beiden Koeffizienten werden wir die Ordnung m zuschreiben. 
AUgemein gesprochen, wird also unsere Intégration bis zur 
Ordnung n einschlieBlich genau sein, welche Werte auch immer 
fur a, a', a" ... a^"^ angenommen werden môgen. Doch steht 
darum nichts im Wege, fur geschickt gewâhlte Werte dieser 
GrôBen die Genauigkeit zu einem hôheren Grad zu erheben. 
So haben wir schon oben gesehen, dafi bei der Méthode von, 
Cotes, d..h. fur 

0/ jL êêÙ II» O 

, a = — , a =— , a =— usw., 
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die Genauigkeit sich von selbst zur Ordnung n + 1 einschliefl- 
lich ausdehnt, so oft n eine gerade Zabi ist. Wenn die Werte 
a, a\ a'\ a'" usw. so gewâhlt worden sind, da6 in der Funk- 
tion T" oder TJ" am Anfang ein Glied fortfàllt oder mehrere, 
leuchtet allgemein ein, dafi die Genauigkeit um ebenso viele 
Stufen liber die Ordnung n hinausgerlickt werden wird, als 
Glieder ausgefallen sein werden. Hieraus wird leicht gefolgert, 
dafi man, da die Menge der Grôfien^ die zur Auswahl stehen, 
w + 1 ist, durch geeignete Bestimmung derselben die Genauig- 
keit immer auf die Ordnung 2 n + 1 eiuschlieBlicb erheben 
kann, und so wird sicb mit Hilfe von n + 1 Gliedem derselbe 
Genauigkeitsgrad erreichen lassen, zu dessen Erlangung es 
nôtig wâre, 2 w + 1 oder 2n + 2 Glieder zu benutzen, wenn 
wir der Méthode von Cotes folgten. 

16. 

Dièse ganze Angelegenheit dreht sich darum, fUr einen be- 
liebigen gegebenen Wert von n eine Funktion T zu ermitteln 
von der Form 

<"+! + a P + a' P-^ + a" r-2 usw. 

derart, daB in dem entwickelten Produkt 

T{t-' + ^t-^+^t-^ + ^t-'+ usw.) 

die Potenzen 

t-\ t-\ <-» ... ^-(«+1) 

aile den Koeffizienten bekommen; oder, wenn man lieber will, 
eine Funktion TJ von der Form 

w«+i + /9 w« + /S' w»»-! + /S" M»-2 + usw., 

deren Produkt mit 

Oui 

von den Potenzen 

w-S u-^, u-^, u"^ ... tt-^'^+i) 
frei wird. Die letztere Art wird etwas einfacher sein: da sich 
nàmlich leicht ersehen lâBt, daB die Koeffizienten von Î7, um 
der vorgeschriebenen Bedingung zu gentigen, intermittierend ver- 
schwinden miissen oder /S = 0, /S" = 0, /S"" = usw. zu setzen 

4* 
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ist, kann fast die Hâlfte der Arbeit schon als erledigt gelten. 
Wir wollen einige einfache Fâlle vorfûhren. 

I. Fur w = mu6 nur der Koefi&zient von <~^ in dem 
Produkt 

verschwinden. Da er = - + a wîrd, haben wir 
£^ = — i- oder T = < — ~ . 

Ebenso folgt Z7=w. 

II. Fur n =: 1 hângt die Bestimmung des T von den beiden 
Gleichongen 

o = y + i-« + «'. 
o = ^ + |« + i«' 

ab, woraus wir ableiten 

a = -l, a'^ + ~ oder T =/«-/ + 4-- 

D 6 

Die Bestimmung der Fnnktion U bringt nur eine einzige Glei- 
chung 



0=| + /î' 



mit sich. Daher 



§'^-\ oder cr=wa^i. 

IIL Fur » =3 2 wird die Funktion T mit Hilfe von drei 
Gleicbungen bestimmt: 






woraus wir gewinnen 



3 

a = — • 

also 



« = -¥' «-5' « =-2Ô' 



^"^ -2^ +5^ 20 
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Zur Bestimmung von TJ geniigt eine einzige Gleichung: 
Daher 



= 1 + 1/5'. 



/S' = — — oder Ï7==w'— — w. 

Doch wollen wir dièses Yerfahren^ welches bestandig lâsti- 
gere Rechnungen herbeifûhrt, hier nicht weiter verfolgen, son- 
dem zur echten Quelle der allgemeinen Lôsung fortschreiten. 

17- 

Wenn ein Kettenbruch vorgelegt ist 

V 



V' 
w -\ 



V" 

w' + — 



V 
w" + 



w'" + usw. * 



so steht fest, daB man fortlaufend bessere Nâherungsbruclie 
durch folgenden Algorithmus findet. Es môgen zwei Beihen 
von GrôBen gebildet werden, F, F', F", F'" usw., TF, W\ PF", 
TF'" usw. durch nachstehende Formeln: 



F =0, 






TF =1, 


7' ^v. 






W ^wW, 


V" =w'F' + 


v'V 




W" ^w'W +v'W, 


F'" =m;"F" + 


v" F', 




FF'" =tt;"TF" +v"TF; 


F"" = w'" V" + 


t)'" F" 




TF"" = ^^;'"TF'"+^;'"TF' 


. Dann wird 


W =0' 








F' V 


> 






F" r 








FF" ~ 

w 


+ 






F'" t; 

IF'" 

to 


-1- 


r' 




^.^ 



usf. 
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AuBerdem steht fest oder wird aus den vorangehenden Glei- 
chungen leicht bestâtigt, dafi 

VW - F' TT = - v, 
usw. ist. Darans ist ersichtlich, daô von der Reihe 



V V V V 



i TXT'lWI't» 1X7"' iA/"" "T" U.8W. 



WW W'W" ^ W"W'" W" W 



y> 

das erste Glied = -^^ ist, 



die Summe der zwei ersten Glieder = 
die Summe der drei ersten Glieder = 
die Summe der vier ersten Glieder = 



W 

V" 



W" ' 

y II' 

yiiii 

us£ Daher wird die Reihe selbst, entweder ins Unendliche 

oder bis zum Abbruch fortgesetzt, den Kettenbruch g) aus- 

drticken. Zugleich erhâlt man hieraus die Differenz zwischen cp 

yi y II y m 
und den einzelnen Nâherungsbrûchen -^^, W^ ' W^ ^^^* 

Aus Formel 33 in Art. 14 der Disquisitiones géné- 
rales circa seriem infinitam* erhalten wir leicht dureh 

Verwandlung von ^ in — die Umformung der Reihe 

tp = vT^ + "q-^"^ + y^""^ + Y^~' + ^8w. 
in folgenden Kettenbruch 



1 
u 

3 



2-2 

3-5 



3-3 

5.7 
u 



4»4 



u — usw. 
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80 dafi man hat 

16 






63 ^^^'^ 



w =^ w = w 



ffieraus bekommen wir fur F, F', F", F'" usw. IF, PF', ÏF", 
W" usw. folgende Werte*: 

F =0, 
V =1, 

F" =M, 

F'" = M* - — , 

F"" =M3_il„^ 

39 ^ 715 15015 ' 

W = l, 

W" = M, 

5 

11 ^ 11 231 ' 
'^ 13 ^ 143 429 

Bei leichter Aufmerksamkeit ist ersichtlich, dafi die ein- 
zelnen F, V\ F", 7'" usw., W, W\ W'\ W"' usw., ganze 
Funktionen der Unbestimmten u bilden, in F^*») w«-^ das 
hôchste Glied wird und die Potenzen u"^-^, u"^-*, w"»-6 ^^g^ 



Digitized by VjOOQIC 



56 Oœuss: Nette Méthode zur nâherungsweisen Intégration 



fehlen, in TT^*») jedoch u"^ das hôchste Glied wird und die Po* 
tenzen u'^'\ w'»-^, w'»~^ usw. fehlen. 

Nach dem oben Bewiesenen aber wird 

^ "■ WW^ "^ ZW'W" "** 3.3.5 W"W"' "^ 3.3.5-5.7 W" W"'" 
2>2.3'3.4>4 

"*" 3.3.5.5.7.7.9 îr""TF'"" "*■ ^^^' 
sein und demnach allgemein 

F^"*^ 2.2.3*3 . . . m*m 

9^~ Pï^(«) ~ 3.3.5.5 ... (2w- l)(2m+ 1) ïr^'"^TF^""*'^> 
, 2.2.3.3 . . . (m + l)(m + 1) 



3.3.5.5 ... (2w+ l)(2w + 3) 17^'"+^) TF ("*•*■ ^^ 
+ USW. 

y(»») 

Wenn man daher œ tt in eine fallende Reihe verwandelt, 

so wird das erste Glied derselben 

2.2.3.3 . . . m-mî^"^^'*'*"^^ 



3.3.5.5 .. . (2m -l)(2m + 1) 
sein. Das Produkt cp TT^*»^ aber wird zusammengesetzt sein ans 
der ganzen Funktion F^"*) und einer unendlichen Reihe, deren 
erstes Glied 

2.2.3.3 . .. mmu~^'^'^^^ 

— 3.3.5.5 ... (2 w - l)(2w 4- 1) ' 

Hiermit ist also von selbst die Funktion U von der 
Ordnung n + 1 gefunden, welche der im vorhergehenden Artikel 
aufgestellten Bedingung genûgt, dafi nâmlich das Produkt (pU 
von den Potenzen 

frei wird. Sie ist nàmlich keine andere als TF^»»+^), und zugleich 
leuchtet ein, daB V gleich F^'^ + i) wird und das erste Glied 
von U" 

^ 2.2.3.3 ... (n+ l)(/^+ 1) .^-(n4.2^ 
3.3-5.5 . . . (2n + l)(2n + 3) 

Wenn demnach fiir 6, h\ V ... 5(") die Wurzeln der Gleichung 
^(n+i) ^ angenommen und die Werte der Koeffizienten i^ 
R\ U" ... -R^*») nach den oben mitgeteilten Vorschriften be- 
stimmt werden, so wird sich unsere Integralformel einer bis zur 
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Ordnung 2n+ 1 aufsteigenden Genauigkeit erfreuen und ihre 
Korrektion wird sehr nahe durch 

1 2.2.3.3 ... (n+l)(n + \) .^^2^^2) " 



22n+2 3.3.5.5 , ,. (2n + l)i2n + 3) 

^ 1.1»2.2.3»3 ... (n-\- l)(n + 1) j^(2n + 2) 

2.6.6.10.10.14 . . . (4n + 2)(4n + 6) 

ausgedrûckt werden. 

18. 

Die Erôrterungen des vorigen Artikels lehren zwar ge- 
eignete Funktionen U fur einzelne Werte der Zahl n finden, 
doch nur nacheinander, indem von kleineren Werten zu grôBeren 
iiberzugehen ist. Leicht bemerken wir aber, da6 dièse Funk- 
tionen sich allgemein ausdriicken lassen durch 

un^i ^ (^ + ^)n n-1 , (n + l)n(n-l)(n-2 ) 3 

2.(2n+l) ^ 2'4{2n + l)(2n-l) 

(n + l)n(n - l)(n - 2)(n - 3)(n - 4) .5 
2.4-6.(2^+ l)(2w- l)t2ti- 3) ^ 

oder auch, wenn wir uns der Bezeichnung F nach der Richt- 
schnur der oben zitierten Abhandlung bedienen, durch 

un^iF(^^^n, -l(n + l), -(n+4-), t*-^) . 

Dièse Induktion wird leicht in einen strengen Beweis verwandelt 
durch eine allgemein bekannte Méthode oder, wenn man es so 
ansehen will, mit Hilfe der Formel 19 in der zitierten Ab- 
handlung.® Die Funktion U kann auch, wenn man es vor- 
zieht, mit umgekehrter Anordnung der Glieder ausgedrûckt 
werden durch 

, 3.5.7 ... (^ + 1) ^/ 1 1 / , ON 3 2\ 

± (n-|-3)(^ + 5)...(2n-fl) '^^i-2^> 2(^ + 3)' 2' ^J 

fiir gerades n, indem das obère oder untere Zeichen gilt, je 
nachdem — n gerade oder ungerade ist, oder durch 



2 
1.3.5 



(n + 2) (^ + 4) . . . (2 n + 1) 



F(-l{n+l). ln+1, 1. «») 



fur ungerades n, indem das obère oder untere Zeichen gilt, je 
2" 



nachdem — (n + 1) gerade oder ungerade ist 
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fachere bevorzugt. Wie dièse sich an die Entwicklung der 
Reihe 

U"^ + — W""^+ -rW^ + U8W. 
6 Ô 

in einen Kettenbruch angelehnt hat, so hâtten wir durch âhn- 
liche Schltisse ans der Entwicklung der Reihe 

in den Kettenbruch 

1 



t- 



1 - 



t- 



_2 



1 - 



2 

10 



3 



1 — usw. 



einen Algorithmus zur Bestimmung der Funktion T fur die 
sukzessiven Werte der Zahl n ableiten kônnen. Zu demselben 
Ergebnis gelangen wir aber durch die Erwâgung, daB T nichts 

anderes ist als — — ^ oder -r-, wenn fiir u geschrieben 

wird 2t — 1, und so wird man die nacheînander fur T anzu- 
nehmenden Funktionen durch folgenden Algorithmus erhalten: 
TT = 1, 

^ 5 ^ 20 ' 
10^ y 2J 8 ^ TÔ:14 4 ^ -«-^« 

+ 7 ^ 7 ^ + 70 

USW. 
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Durch Induktion ergibt sich hieraus allgemein 

1.(2^ + 2) ^ 1.2.(2n+ 2)(2n+ 1) 



1.2.3.(2 n + 2)(2fi +l).2n 

oder 

T=/«+iir(-(n+l), ~(n+l), -2(n>l), 0> 

und es ist leicht, dieser Induktion die Eraft eines Beweises zu 
verschaffen. Wenn es angenehmer ist, kann T auch mit um- 
gekehrter Anordnnng der Glieder durch 

± .r.-;..v...',:::v ^('-^-^' -'"+■'-''') 

ausgedrtlckt werden, wo das obère Vorzeichen fttr ungerades n 
gilt, das untere ftLr gerades. Auf âhnliche Weise wird endlich 
allgemein gefunden T" gleich dem Produkt aus 

1.1.2.2.3.3 ... (7^ + l).(n + l) _f„^2) 

2.6.6.10-10.14 . .. (4n + 2).(4n + 6) 

und der unendlichen Beihe 

- (n + 2)« 1 (n + 2)«(n + 8)« , 

"*" l.(2n+4) "^ 1.2.(2n + 4)(2n + 5) 

(n + 2)«(n + 3)«(n + 4)« ^.3 _^ ^^^^ 



' 1.2.3.(2fi + 4)(2n + 5)(2n + 6) 
oder 

F{n + 2, n + 2, 2 n + 4, /"i). 

30. 

Da in der Funktion U die Potenzen w»», tt"-^^ ^n-4 ^g^^ 
fehlen, so werden von den Wurzeln der Gleichung 17 = stets 
je zwei entgegengesetzt gleich sein. Dazu muB man fÛr einen 
geraden Wert von n noch die singulâre Wurzel gesellen. 
Nach Auffindung der Wurzeln wird man die Werte der Koeffi- 
zienten R, R\ R" usw. gemâB der Méthode des Artikels 11 er- 
halten durch eine ganze Funktion von u, die flir einen un- 
geraden Wert des n von der Form 

y t^n-l ^ y'î^«-3 ^ j/"wn-5 ^ ug^, 

sein wird, fUr einen geraden Wert aber, wenn der der Wurzel 
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w = entsprechende Koeffizient ausgeschlossen wird, von der 
Form 

/W**-2+ J^'w**-^ + /"W*»-^ + U8W. 

Das Beispiel des Ârtikels 12 zeigt eben dièse Reduktion fiir 
w = 6. Offenbar entsprechen also entgegengesetzten Werten 
von u stets gleiche Koeffizienten. Ubrigens kann in dem Falle, 
wo n gerade ist, der der Wurzel w = entsprechende Wert 
leicht allgemein a priori bestimmt werden. Man gewinnt diesen 

Koeffizienten, wenn man in j^jjT einsetzt w = 0. Den Wert 

des Zâhlers U' fiir w = haben wir bereits in Artikel 18 mit- 
geteilt, der Wert des Nenners aber wird ebendaher 

= , 3'5.7 ... (n-f- 1) ^ 3' 3«5-5'7 »7 . . . (n + 1) (n -f 1) 

""^ (w+ 'i){n + 5) . . . i2/ï + 1) ■"* " 8.5.7.9«11 . . . (2n + Ij 

sein und voUends der gesuchte Koeffizient 



"" \'3^.7«9 . .. (w + l)j • 



21. 

Die ganze Funktion von u, welche die Koeffizienten R, 
R\ R" usw. darstellt, kànn in dem hier behandelten Falle auch 
unabhàngig von der allgemeinen Méthode des Artikels 11 
folgendermaBen ermittelt werden. Durch DifiFerentiation der 
Gleichung 

dann durch die Substitution 

dq> _ 1 
du "~ 1 — w* 

und Multiplikation mit U^{u^^l) erhalten wir 

Die Glieder dieser Gleichung zur Linken bilden ofifenbar eine 
ganze Funktion von u: deshalb mtissen in dem Teile zur Rech- 
ten die Koeffizienten von u mit negativen Exponenten not- 
wendig sich selbst zerstôren. 
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ÎT 
Doch — — bringt eine unendliche Reihe hervor, die mit 



du 
dem Gliede ' 

_ / 1.2.3.4 .,.(n + l) \^^-r2n^4^ 
\1.3.5.7 ... (2/^ + l)/ 

beginnt, und wenn man dièse also mit [u^ — 1) U^ multipliziert, 
so wird nichts anderes herauskommen konnen als die konstante 
GrôBe 

_ / l«2-3»4 ... (n + 1) \^ 
\ 1.3-5.7 . . . (2n+ 1) / * 

Hieraus schlieBen wir*, dafi 

(^ - ^) ^ ^iT + i 1.3.5.7. ..(2n+l)j 

durch U teilbar ist. Deswegen wird mit der gebrochenen Funk- 

tion i^Tj\ y welche die Koeffizienten J% R\ i2" usw. liefert, 

\dû] 
die ganze Funktion 

__ / 1.3.5.7... (2n-H) ^.V . ., _ 1) 

V 1.2.3.4 ... (n 4- 1) ; ^ ^ 

gleichwertig sein. Statt dieser Funktion, die von der Ordnung 
2 n + 2 ist und oflFenbar bloB gerade Potenzen von u enthàlt, 
-wird der Rest aus ihrer Division durch U genommen werden 
kônnen, welchem die Ordnung n oder » — 1 zukommt, je nach- 
dem n gerade oder ungerade ist, Wenn wir aber in dem er- 
steren Falle lieber den Koeffizienten, der der ,Wurzel w = 
entspricht, ausschlieBen woUen, so werden wir statt jener Funk- 
tion deren Rest aus der Division durch — wâhlen, welcher sich 
nur ZUT Ordnung n — 2 erheben wird. 

23. 

Damit zur Hand sei, was zur Anwendung der bisher aus- 
einandergesetzten Méthode erfordert wird, erschien es gut, fQr 



* „Zugleich wird hieraus der Beweis entnommen, daB U nicht einen 
stimmten gemeinsamen Teiler mit — — ha 
die Gleichung Î7 = keine gleichen Wurzehi." 



iinbestimmten gemeinsamen Teiler mit — — haben kann und demnach 
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die sukzessiven Werte der Zahl n die numerischen Betràge so- 
wohl der GrôBen a, a', a" usw., als auch der Koeffizienten R, 
E', R'' usw.^ auf 16 Stellen berechnet^ zusammen mit deren 
Logarithmen bis zu zehn Stellen beizafûgen. 

I. Ein Glied, n = 0. 

a = 0,5, 
R^ 1. 

Die Korrektion der Integralformel sehr nahe = -r^ U\ 
IL Zwei Glieder, n=l. 

T = <»-< + -!-, r'=<_i-, 

a =0,211324865i051871, 
a' = 0,7886751345948129, 

Die Korrektion sehr nahe = —^^r L"'\ 

III. Drei Glieder, n = 2. 

J - « - 2 ' + 5 < - 20 ' ■' - ' * + 60 

a =0,1127016653792583, 

a' = 0,5, 

a"= 0,8872983346207417, 

Die Korrektion sehr nahe = -rr— r- L^'. 

AOUU 
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IV. Vier Glieder, n = 3. 

T =.t*-2t^ + \t'-\t + ^, 

J -< -yî + 21 ' 84 ' 

a =0,0694318442029754, 

a' =0,3300094782075677, 

a" «0,6699905217924323, 

a'" = 0,9305681557970246, 
R =JÎ"' = 0,1739274225687284, log = 9,2403680612, 
R'^R" =0,3260725774312716, log = 9,5133142764. 

Von diesen Koeffizienten der allgemeine Auadruck: 

35 i , 11 
-W** +18-' 

die Korrektion sehr nahe = -^^ I^^"^- 

V. Fanf Glieder, n = 4. 

T=t*-2t + 36 ' -36»+ 75600' 
a = 0,0469100770306680, 
a' =0,2307653449471585, 
a" =0,5, 

a'" = 0,7692346550528415, 
a"" = 0,9530899229693320, 

KowalewBki, Newton usw. 
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R =5""= 0,1184634425280945, log = 9,0735843490, 
R' =,R"' = 0,2893143352496832 9,3789687142, 

R"=P-= 0,2844444444444444 9,4539974559. 

225 

Der allgemeine Aasdruck dieser Koeffizienten, unter Aus- 
schluB von R'': 

400 ^ "*" 3600 ' 

die Korrektion sehr nahe = ^kzttt L^- 

O0Od44 



VL Sechs Glieder, n = 5. 



5_ 
231 



m' ^ ,5 _ i.4_ ''l<S_i9,2 , ^, 1_ 

2 33 22 ''' :J20 1320 ' 

a =0,0337652428984240, 

a' =0,1693953067668678, 

a" =0,3806904069584015, 

a'" =0,6193095930415985, 

a"" =0,8306046932331322, 

a'"" = 0,9662347571015760, 
B =iî""'= 0,0856622461895852, log = 8,9327894580, 
R' =R"" = 0,1803807865240693 9,2561902763, 

R" = R"' =0,2339569672863455 9,3691359831. 

Der allgemeine Âusdmck der Koeffizienten: 

800 75 ^ 96 ' 

die Korrektion sehr nahe = L^^ . 
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VIL Sieben Glieder, n = 6. 

TT 7 21 5 , 105 o 35 

13 ^ 143 429 ' 

TTf A 50 4 , 283 2 256 

2 ^ 13 52 ^ 143 286 ^429 3432' 

^156 78 ^2860 4290 ^ 240240 ' 

a = 0,025 446 043 828 620 2 , 

a' =0,1292344072003028, 

a" =0,2970774243113015, 

a'" =0,5, 

a"" = 0,7029225756886985, 

a'"" =0,8707655927996972, 

a""" = 0,9745539561713798, 
B = R — = 0,0647424830844348, log = 8,8111893529, 
R' =R""' =0,1398526957446384 9,1456708421, 

R" = R"" =0,1909150252525595 9,2808401093, 

fi"'= .^ = 0,2089795918367347 9,3201038766. 

Der allgemeine Aasdrnck dieser Eoeffizienten, unter Ans- 
scbloB von R": 

1859 . 1573 , , 7947 



16800 29400 ~ 39200 ' 

die Korrektion sehr nahe = 175 679 360 ^'^* 



33. 

Zum SchluB voUen wir die Leistungsf&bigkeit unserer Mé- 
thode Yor Augen stellen durch Berechnong des Integralwertes 

dx 



/i 



log a; 
von X = 100000 bis x = 200000. 
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I. Au8 einem Gliede erhalten wir 



II. Au8 zwei Gliedern ergibt sich 



III. Au8 drei Gliedern ergibt sich 



IV. Ans vier Gliedern ergibt sich 



JR A 


= 8390,394608 


AR A 


= 4271,810097 


AR' A' 


= 4134,144502 



V. Au8 fûnf Gliedern ergibt sich 



VI. Aus sechs Gliedern ergibt sich 



Sa. = 8405,954599 
AR A = 2390,572772 
AR' A' =3729,064270 
AB" A" =2286,599733 



Sa. =8406,236775 
AR A = 1501,957053 
AR' A' = 2763,769240 
AR" A" = 2711,454637 
AR'" A'" = 1429,062040 





Sa. 


= 8406,242970 


AR 


A 


= 1024,879445 


AR' 


A' 


= 2041,833335 


AR" 


A" 


= 2386,601 133 


AR'" 


A'" 


= 1980,509616 


AR""A"" 


= 972,419588 




Sa. 


= 8406,243117 


AR 


A 


= 741,912854 


AR' 


A' 


= 1545,757256 


AR" 


A" 


= 1976,737668 


AR" 


A'" 


= 1950,466223 


AR"" A"" 


= 1488,588550 


AR^ 


^^ 


= 702,780570 



VII. Aus sieben Gliedern ergibt sich 



Sa. = 8406,243 121 
AR A = 561,1213804 
AB' A' = 1202,0551998 
AR" A" = 1621,6290819 
AR'" A'"= 1753,4212406 
AB"" A""= 1584,9790252 
JJ?v A^ = 1152,0681116 
J^vi^vi= 530,9690816 



Sa. =8406,2431211 

Nach den Berechnungen des berûhmten Bessel ist der Wert 
desselben Intégrais = 8406,24312 gefunden worden. 



Digitized byCjOOQlC 



Jacobi: Ûber Oauss' neue Méthode usw. 69 



Uber Gauss' neue Méthode, 
die Werte der Intégrale nâherungsweise zu finden. 

(Von Cari Gustav Jacob Jacobi, 1826.) 



1. 

In den Principiis von Newton^ liest man eine Méthode, 
wie man durcb eine Anzahl gegebener Punkte eine parabolische 
Kurve legen kônne. Dièse Aufgabe erscheint analytisch als 
Interpolationsproblem, ans mehreren Gliedern einer Reihe das 
allgemeine zn finden. Es ist der bekanntere Fall, wenn die 
Intervalle der Ordinaten der gegebenen Punkte gleich groB 
sind, oder analytisch ausgedruckt, wenn die Werte des reihenden 
Eléments, fur welche auch die Werte der entsprechenden Glie- 
der der Reihe gegeben sind, eine arithmetische Progression 
bilden. Aber der élégante, mit Unrecht weniger gekannte 
Algorithmus, den Newton gibt, erstreckt sich schon auf den 
allgemeineren Fall, wenn jene Intervalle der Ordinaten der ge- 
gebenen Punkte oder jene Werte des reihenden Eléments irgend 
beliebige sind. Newton hat hiervon eine Anwendung auf die 
Quadraturen gemacht. Durch mehrere Punkte der zu quadrie- 
renden Kurve, fur welche die Ordinaten berechnet worden sind, 
legt er die parabolische Kurve, und deren Quadratur zwischen 
denselben Grenzen, zwischen denen die gegebene Kurve qua- 
driert werden soUte, gibt einen Nâherungswert. 

Newton hat von jenem Interpolationsproblem und seiner 
Anwendung auf die Quadraturen ferner in einem Traktâtchen 
gehandelt, welches Methodus diflferentialis betitelt ist und zuerst 
in der Amsterdamer* Ausgabe seiner Principia, vom Jahre 1723, 
nebst anderen Abhandlungen angehângt gefiinden wird. Hier 

* Von dieser Ausgabe ist die Kuriositât zu erzâhlen, daB aie auf 
Kosten des beriilimten Philologen Richard Bentley veranstaltet worden 
ist, der in seinen englischen und lateinischen Predigten oft die Principia 
seines genauen Freundes Newton anpries, als ein Bollwerk gegen die 
Irreligiositât und eine OfiFenbarung der GrôBe Gottes. 

(Anmerkung Jacobi s.) 
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rat er unter anderem zum Behuf der leichteren Berechnung der 
Intégrale ftir jede Zahl der berechneten Ordinaten, deren Inter- 
valle er gleich gro6 annimmt, Tafeln anzufertigen, von denen 
er auch selbst einen Anfang gibt, welchen hemach Roger 
Cotes in seiner Harmonia mensurarum fortgesetzt hat 

Aber Gauss hat in den Gôttinger Kommentarien gezeigt, 
daB man durch schickliche Wahl der Abszissen, fiir welche die 
Ordinaten berechnet werden, den Grad der Nâherung auf das 
Doppelte treiben kann; und da solche Bestimmung unabhângîg 
von der Natur der zu quadrierenden Kurve geschieht, so ist es 
moglich, auch nach der so yervoUkommneten Méthode Tafeln zu 
verfertigen, von denen auch Gauss eine Probe gegeben hat. 
Gauss gelangt zu seinen Resultaten auf dem Wege einer 
schwierigen Induktion, die durch die sogenannte Kàstnersche 
Méthode, wenn etwas fur die Zahl n gilt, es auch fur die Zahl 
n + 1 zn erweisen, zur AUgemeinheit erhoben werden kann. Es 
ist also noch ein direkter Beweis zu wiinschen. Die groBe Ein- 
fachheit und Eleganz der Gaussischen Eesultate làBt einen 
einfachen Weg vermuten. Auf einem solchen einfachen und 
direkten Wege zu jenen Resultaten zu gelangen, mit denen 
Gauss die Wissenschaft bereichert hat, ist der Zweck dieser 
Abhandlung. 



Es sei das Intégral i ydx zwischen den Grenzen x = und 

X = 1 zu nehmen. Andere Grenzen werden leicht auf dièse 
zurûckgefûhrt. Es seien ferner die Werte von x, fiir welche y 
bekannt ist, 

ce', a' y a"\ . . . , a^"^ 

so daB, wenn man y = /"(x) setzt, die entsprechenden Werte 
von y werden: 

/•(«'), /•(«"), /•(«'"), ..., A«^"')- 

Man bilde das Produkt 

[x - a)[x - a")[x - «'") ... [x- a'^'') 

und nenne es (p{x)y so hat man, wenn y = f[x) eine ganze ratio- 
nale Funktion vom {n — l)ten Grade ist, durch Zerfâllung in 
Partialbrûche^: 
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A?) _ /(«') ■ ((»") , fi"'") 



wo wir mit ç? '(«("*)) den Wert von 

fur rr = «(*"> bezeichnen* Vermittelst dieser Formel findet man 
durch Multiplikation mit 97 (rr) sogleich y aus den spezielien 
Werten fiir 

x=^ a\ x= a" y X ^ a'", . . . , a? = a^^'K 

Ubersteigt aber y den [n — l)ten Grad, so gibt der Ausdruck 
zur rechten Seite des Gleichheitszeichens, welchen wir O nennen 

wollen, nur den echten Bruch, der in dem unechten -^-f^ steckt, 

so daB, wenn f{x) z. B. vom {n +^)ten Grade ist, und man 

f{x)^ U+ V,(p[x) 

hat, wo U hôchstens vom [n — l)ten, V vom ^ten Grade ist, 

Entwickelt man O und den Bruch -^ nach den absteigenden 
Potenzen von x, so enthâlt 

G= ^ 



q)(x) 

die negativen, V die positiven Potenzen von ar, die sich in der 
Entwicklung von ' ^ befinden. Setzt man daher 

f[x) = a + ax + a"x^ + . . . + a<»^x^» + a<^+^^ir'»+' 

und 

1 A' , A" , A'" , , ^("+^> , 

^ = -^ + ^^ + ^^+-- + "-;^+ usw., 

SO ûndet man 

+ ... + a(2«-i)(^'a;«-i + A" x^-'^ + ... + ^^"0 + ^sw. 
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8. 

Newtons Nâherungsmethode besteht darin: statt y^f{x) 
die Funktion 

U= Q.(p[x) 

zu substituieren. Der Fehler oder die Diflferenz der Intégrale 
der gegebenen und substituierten Funktion wird dann 

J = 1 ydx — I Udx = I y (^) Vdx. 

Es wird jetzt die Aufgabe gestellt, die GrôBen 

a', «", a'", ..., a'^ 

80 ZU bestimmen, da£ der Fehler J môglichst gering oder 
die Nàherung môglichst genau werde. In den Fâllen, wo die 
Nàherungsmethode mit Gltick angewendet werden soU, mtissen 
die Koeffizienten der fur y gesetzten Reihe rasch abnehmen. 
Je mehr daher von den ersten Koeffizienten dieser Reihe, 
welche die hauptsâchlichsten sind, in dem Âusdruck fiir den 
Fehler J verschwinden, desto kleiner wird er im allgemeinen 
und desto grôBer die Nàherung. Da nun schon, was auch die 
GrôBen 

a, a", a'", . . ., a^^ 

waren, im Ausdruck fiir 



J = l (p{x). Vdx, 



me aus dem fur F gefundenen Ausdrucke erhellt, die Koeffizienten 

a, a', a", ..., a^*»-!) 

nicht mehr vorkommen, so wollen wir rermittelst schicklicher 
Bestimmung jener GrôBen auch noch die mit 

behafteten Glieder verschwinden machen, wodurch ein doppelter 
Grad der Nàherung erreicht wird. Es wird dièses immer môg- 
lich sein, da die Zabi der willkiirlichen GrôBen und der zu er- 
ftillenden Bedingungen dieselbe ist. Man sieht sogleich aus dem 
fur V gefundenen Ausdruck, daB hierzu eine solche Bestimmung 
von q>{x) erfordert wird, daB die Intégrale 

j(p{x)dx, j x^{x)dx, jx^(p{x)dx, . . ., j x"^^^ (p{x)dx 
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zwischen den Grenzen rr = und x ^ ly zwischen denen das 
Intégral jydx genommen werden soU, yerschwinden. Dièse 
Bestimmung ist jetzt die Aufgabe. 

4. 

Es lâfit sich darch eine bekannte Redukidonsformel das 
Intégral 

i x^(p {x) dx 

auf die vielfachen Intégrale von cp{x) znrtickfiihren. Man bat 
nâmlich allgemein: 

j uvdx = u l vdx — jrfw/ vdx, 
l du 1 vdx ^ du l vdx — l d^u j vdx, 
l d^u I vdx -= d^u \ vdx — \ d^u I vdx, 

d^u j vdx = d^u l vdx — l d^'^^u l vdx, 

wo man jede Formel ans der vorhergehenden erhâlt, indem man 
-^ statt u und (vdx statt v setzt» Hieraus folgt sogleicb: 

l uvdx ^ u j vdx — du l vdx+d^ut vdx 

/m -|-1 ^% y^TO +1 

vdx + (- Ij'^+i / d'^+^u j vdx. 

Setzt man 

w = «*" , V = (p{x), 

so erbâit man bieraus: 

j x^^{x) dx = a?"* I gp {x) dx — mx^~^ I cp (a?) dx^ + 

+ m (m — 1) aj*" - 2 i ç^^i^x) dx^ — . . . 

(- l)"*m(m ~ \)[m - 2) ... 1 T"* 9?(a:)da;™+^ 
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Gibt man dem m nacheinander die Werte 

0, 1, 2, 3, ..., n- 1, 
so erhalt man: 

l q>{x)dx = j ^ (^) ^^ y 

lxq){x)dx ^ X l <p{x)dx — l <p{x)dx^, 

lx^(p{x)dx =^ x^ I <p{x)dx ^ 2x l fp{x)dx^ + 2j (p{x)dx^, 

j x'*-^(p{x)dx = x""-^ j (p{x)dx — (w — l)a;"-2 l (p{x)dx^ 
+ (w-l)(w- 2)a;"'-3 / (p[x)dx^—.., 

{-l)n-i{n^l){n -2).. . 1 r(p{x)dx''. 

Dièse Formeln sind bekannt Man sieht aus ihnen, dafi^ 
wenn 

I (p {x) dx, l x(p{x)dx, I x^(p{x)dXf . . ., i x"~^ ^ (a?) dx 

zwischen gewissen Grenzen verschwinden sollen, zwischen den- 
selben Grenzen auch 

(p[x)dXj I (p{x)dx^, l cp{x)dx^y . • •> / (f{x)dx'' 
verschwinden mtissen und umgekehrt. ^ 

5. 

Unsere Aufgabe ist also jetzt darauf zuriickgefùhrt, die 
Funktion y (a;) so zu bestimmen, daB ihr Ites, 2 tes, 3 tes, . . ., 
wtes Intégral zwischen den Grenzen a? = und a; = 1 ver- 
schwinden; d. h., wenn man die aufeinanderfolgenden Intégrale 
bis zum nten so bestimmt^ daB sie fur a; = verschwinden, so 
sollen sie auch fur a; = 1 verschwinden. 

Man setze 

I (p{x)dx^ = n{x)t , 
die aufeinander folgenden Intégrale so bestimmt, daB jedes fiir 
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x == yerschwindet, so kann man jetzt die Aufgabe so aus- 
drûcken, eine Funktion n [x) zu finden, die flir rr = und fur 
a; = 1 zugleich mit ihrem Iten, 2ten, 3ten, ..., (w — l)teii 
Differentiale yerschwindet. Dièses erheischt, dafi die Funktion 
n[3^ die Faktoren a;" und (x — 1)" habe, und umgekehrt, jede 
Funktion, die den Faktor x" {x — If hat, erfûUt die verlangten 
Bedingungen. Es mu6 daher gesetzt werden: 

n(x)r^x''[x- ITM. 
Da nun 

(p{x) = (a; - a'){x - a"){x - a'") . . . (x - a^% 

also eine ganze rationale Funktion von der wten Ordnung ist, 
so ist 



(p{x)dx^ 



eine ganze rationale Funktion von der 2wten Ordnung, woraus 
folgt, daB M fur unsern Fall eine Konstante ist. Auf dièse 
Weise erhàlt man* 



rfa;" 2n ' 1-2-2 n {2 n - 1) 

n«(n~l)«(n-2)^ 3 .„ _n(n-lKn - 2) .^ 

l'2-3'2n(2n-l)(2n-2) "T-'-V ) 2n{2n-'l)(2n-2)...(n-h'i)^ 

wo 

M = 



2 w (2 » - 1) (2 n - 2). . .(^ + 1) 

gesetzt worden ist. 

Die Wurzeln der Gleichung y(a;) = 0, ftir (p{x) den eben 
gefundenen Ausdruck gesetzt, geben dann die GrôBen 

a', a", a'". . . ., «<«>, 

so bestimmt, dafi der Grad der Nâherung der môglichst grôfite 
sei. Da aus der Lehre von den Gleichungen bekannt ist, daB, 
wenn die Wurzeln einer Gleichung n{x) = aile reell sind, 
auch aile Wurzeln der Gleichung 

d'Uni x) ^ Q 
dx"^ 
reell sind und zwischen den Wurzeln jener Gleichung liegen, 
so folgt hieraus, da die Wurzeln der Gleichung n(x) = oder 
der Gleichung 

a;"(a;- 1)~ = 
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aile reell sind^ und zwar n von ihnen = 0, die anderen = 1, 
da6 auch die Wurzeln der Gleichung y (a;) = oder die GrôBen 

a\ a'\ a'\ . . ., a<"» 

aile reell sind und zwiscben und 1 liegen^ wie es auch G au s s 
in den berechneten Beispielen gefunden bat. 



In unserer (§ 4) gefundenen Formel: 
j uvdx = u l vdx ^ du j vdx + d^u l vdx 

seize man 

m=7W — 1, w=F, V = q> {x), 

80 erbëlt man, da die n ersten Intégrale von v = g){x) zwiscben 
den Grenzen a; = und x = 1 rerscbwinden, und 



(p [x) dx"" — ^ 



2n(2n - 1) . . . (fi + 1) ' 



rfa?"- 



> 



welcbes Intégral zwiscben den Grenzen ce = und x ^ 1 
zu nebmen ist 

Man seize femer in der angefubrien Formel w = /'»+i und 
es verscbwinde t flir x ^ l, so wird aucb 

du d^u d^u d'"^u 

fiir X = l verscbwinden. Es seien femer die Intégrale 

jvdx, I vdx^, I vdx^, ..., 1 vdx"^-^^ 

so genommen, da6 sie insgesamt fur x = verscbwinden; so 
verscbwinden 

ujvdx, du j vdx, d^u I vdx, . . ,, d^^ul vdx 

zwiscben den Grenzen a; = und x = L Man erbâlt demnach 
zwiscben den Grenzen a; = und a; = /: 

Çuvdx = ff^-^^vdx = (— 1)™+^ fd'^+i ^™+i r vdx. 
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Setzt man jetzt 

d^ V 
t = l—x, Z = l, m = w — 1, V = x'^ , , 

80 erhâlt man zwischen den Grenzen a; = und a; = 1 : 



1.2.3... wl"^ x'^d'^Vdx. 



Man erhâlt auf dièse Weise: 

21 = : I x^d^Vdx. 

2n(2»- 1) ... (» + 1) J ^ "'»' «*.*', 

wo die aufeinander folgenden Intégrale so zn nehmen sind, dafi 
aie fur a; = verschwinden und, nach beendigter Intégration, 
rr = 1 zu setzen ist Unter dieser Form ist der Fehler J am 
leichtesten za berechnen. 

7. 
Vermôge des § 2 findet man: 

~d^ '^ 

+ at2u+i)((^ + i),^^_l) ...2ul'a;+w(w-^l)(^t-2)... 1^'') 

+ a(2«+2) ((n + 2) (w + 1) w ... 3.i'a;2 + (^ + 1)^ (^ _ i) . . . 2 A' x 

+ w(w-l)(n-2)... 1^'") 

+ a(2n+3. ((^_f.3)(^4.2)(n + l)...4^'x3 + (w + 2)(n + l)w...3^'V 
+ (n+l)7î(n-l)...2^'"a; + w(w-l)(n-.2)...l^iv) 

usw. 
Hieraus ergibt sich 

. ig.2«.3g . . . w^ 

"■ (w + l)2(n + 2)2 ... (2 w)«(2 w + 1) 
a(2«) ^' + a(2«+i) (I^^Ç ^' + ^ j + 
«(2«+2) / (n+l)Mn + 2)2 , ^.^1)» „ ,A 



.(2n+3) / (n +l)^(/i + 2)M7i+3)» ^ _ 

\^1.2.3.(2n + 2)(2w+3)(2w + 4) ^ !« 



2-(2n + 2)(2w+3) 



USW. 
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Dièse ersten Glieder des Fehlers J kônnen zur Eorrektar 
dienen. Die GrôBen A', A", A'", A^^ usw. bilden eine wieder- 
kehrende Reihe, da sie aus der Entwicklang des Brucbs 



(p(x) 

1 



2n 1.2.2w(2n-l) 1.2.3-2 n(2n - l)(2n - 2J 

+ ...(-l)n ^n~l)(n-2) ... l 



2n(2n -l)(2n - 2) . . . (n +1) 

entstanden sind, welche wir (§ 2) 

A' , A" . ^'" , A^^ , 

gesetzt hatten. Sie werden durch die Gleichungen gefunden: 

2n ' 

''"^ 1.2.2n(2n-l) ^ 2n^^ ' 

= ^/ n'{n-\f(n-2)^ ,, n«(n-l)« „, n^ _ i^ 

1.2.3.2n(27i-l)(2n-2) 1.2.2n(2n-l) "*" 2n 

USW. 

Dièse Resultate stimmen genau mit den von G au s s ge- 
fundenen tiberein. 
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I. Newton. 

Die Newtonsche iDterpolationsmethode kommt schon in 
seinem bertihinten Werke Philosophiae naturalis principia mathe- 
matica (1687) vor. Da auf die betreflfende Stelle auch in der 
Abhandlung von Cotes Bezug genommen wird, so sei sie hier 
voUstândig wiedergegeben: 



Drittes Bach. Fûnfter Hillssatz. 

Eine krnmme Linie parabolischer Natur zu finden, 
die durch irgendeine Anzahl gegebener Punkte hin- 
durchgeht. 

Jene Pankte seien A, B, C, D, E, F, ... und m an fàlle 
von ihnen anf eine gegebene Gerade beliebiger Lage die Lote'*' 
AH, BJ, GK, DL, EM, FN, . . . 

Erster FalL Wenn die Intervalle HJ, JE, KL, . . . der 
Punkte H, J, K, L, M, N gleich sind, so bilde man von den 
Loten AB, BJ, GK, ... die ersten Diflferenzen** h^ b^, h^, b^, 
b^, . . . , die zweiten c, c^, e^, c^, . . . , die dritten d, d^, d^, . . . , 
80 daB also 

AH'-BJ^b, BJ- GK=b^, GK-DL^b^y 
DL + EM^b^, -EM+FN^b^, ..., 
femer 6 — Jg = c ist usw. 



* Im Original steht „beliebig viele Lote". 

** Newton lâBt den Index 1 fort. Er schreibt ûberdies die andem 
Indizes voran: 2 b, Sb, 4 b, bb, ... 
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Nun werde irgendein Lot RS errichtet, das eine Ordinate 
der gesuchten Eurve darstellt. Um die Lange desselben zu 
iinden, nehme man an, dafi die Intervalle HJ, JK, KL, 
LMj ... gleich der Einheit seien, and setze: 



AH^ 



\r.{+SL)^s, \s>[+SM)^t. 



So geht es also fort bis zu dem vorletzten Lot ME, und man 
hat dabei die Glieder ES, JS, . . . , die von iS ans nach A zu 



/l\''' 



ô 



^c 



M >v 



/y j s /r 



à 6, 4 4 6^ 

^ ^z ^s ^. 

c/ < y. 



£\ 



Fig. 6. 

liegen, mit negativem Zeichen zu versehen, die Glieder SK, 
SL, . . ., die auf der andern Seite des Punktes S liegen, mit 
positivem Zeichen. Hat man die Zeichen richtig beachtet, so 
vrird 

RS ^ a + bp + cq + dr + es + ft + ... 
sein. 

Zweiter Fall Wenn die Intervalle HJ, JK, ... der Punkte 
H, J, K, L, . . . ungleich sind, so nehme man von den Loten 
AH, BJ, GK, ... die ersten Diflferenzen, dividiert durch die 
Lotintervalle, und nenne sie h, b^, b^, b^, b^, . . ., die zweiten, 
dividiert durch je zwei Intervalle, nenne man c, c,^, Cg, c^, . . ., 
die dritten, dividiert durch je drei Intervalle, d, d^, d^, . . ., die 
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vierten, dividiert durch je vier Intervalle, e, e.^, ... usf. Hier- 
bei ist also 

, ÂH-BJ __ BJ^CK 

EJ ' ^^^ JK ' 

^3 ^ KL ' ^* "" ^^ ' " ' 
ferner 

dann weiter 



^ " HK ' ^^^ JL ' "3 - KM 



J ^ ~ ^2 J ^2 ~ ^8 



EU 






Sind die Differenzen gefunden, so seize man 

AU^a, -ESr=.p, P'{-JS) = q, q'{+SK)^r, 

r.{+SL)^s, S'{+SM)^t, 

So geht es also fort bis zum vorletzten Lot ME. Dann 
wird flir die Ordinate BS die .Gleicbung gelten: 

ES = a + bp + cq + dr + es + ft + ... 

Folg&rung, Hiervon ausgehend kann man die Flâchen .aller 
Kurven beliebig genau ermitteln. Findet man nâmlich von 
einer zu quadrierenden Kurve irgendwieviele Punkte und denkt 
sich durch dieselben eine Parabel gelegt, so wird die Flàche 
dieser Parabel ziemlich genau mit der Flâche jener zu qua- 
drierenden Kurve ubereinstimmen. Die Parabel làBt sich aber 
nach wohlbekannten Methoden geometrisch quadrieren. 



Newton benutzt den obigen Hilfssatz, um „aus irgendeiner 
Anzahl beobachteter Ôrter eines Kometen dessen Ort fiir einen 
beliebig gegebenen Zwischenpunkt*' zu ermitteln (was auch Cotes 
am Anfang seiner Arbeit erwâhnt). Erst 1711 kam er aber 



" *~ THJf ^ "^ ~ "771/' ûaben wir hmzugefugt. 



LM "^" EM 

Kowalewski, Newton usw. 



Digitized by VjOOQIC 



82 AnhaTig: L Newton 



dazu, durch Jones eine selbstândige AbhandluDg liber seine 
Interpolations- und Quadraturmethode heransgeben zu lassen 
unter dem Titel: Methodus differentialis. 

Zwischen diesen beiden Newtonschen Auslassungen liber 
die genannte Méthode besteht aber ein erheblicher Unterschied. 
Die Hauptformel des Lemma V, welches wir oben wieder- 
gegeben haben, kommt in der Methodus differentialis 
nicht vor. Dièse Formel 

ES =^ a + bp + cq + dr + es + ft+ ... 

ist gerade das, was man jetzt gewôhnlich als „Newtonsche 
Interpolationsformel*' bezeichnet. Es wird bei dieser Formel,^ 
kurz gesagt, ein Eand des Schémas der Differenzenquotienten 
als gegeben betrachtet. In der Methodus differentialis dagegen 
ist die Mittellinie des Schémas ausgezeichnet: die auf ihr 
liegenden Differenzenquotienten bzw. die arithmetischen Mittel 
aus zwei angrenzenden Differenzenquotienten werden als bekannt 
angesehen. Aus den Betrachtungen in Propositio II der Me- 
thodus differentialis kann man die Interpolationsformel des 
Lemma V ohne weiteres entnehmen, da in Propositio II gerade 
mit den Differenzenquotienten eines Randes operiert wird. 
Trotzdem lâBt Newton seine altère Interpolationsformel ganz 
unerwâhnt. Es liegen allerdings 24 Jahre zwischen beiden 
Publikationen. Môglicherweise hat Newton bei der Àbfassung 
der Methodus differentialis nicht mehr an die altère Formel 
gedacht. In einer Arbeit von Stirling, die in den Philo- 
sophical Transactions von 1719 unter dem Titel Methodus diffe- 
rentialis Newtoniana illustrata erschien, werden bei de Formeln 
entwickelt (freilich nur fur gleiche Intervalle). Es wird dort 
zugleich mit Nachdruck betont, da6 Newton schon 1676 sein 
Interpolationsverfahren in Vorlesungen behandelte. Vielleicht 
sah es Newton ganz gern, daB auf solche Weise auch seine 
altère Formel wieder ans Licht gebracht wurde. Der Erfolg 
war leider der, da6 man spâter die zweite Newtonsche Formel 
vielfach auf Stirlings Konto setzte. So schwer làBt sich in 
der Wissenschaft Eigentum sichern. 

Moritz Cantor bespricht in seiner monumentalen Ge- 
schichte der Mathematik (Bd. III, 1898, S. 358 ff.) ausfûhrlich 
die Methodus differentialis und tut so, als behandelte sie nur 
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die erste Newtonsche Interpolationsformel, wâhrend, wie wir 
schon bemerkten, gerade nur die zweite darin vorkommt. 

Der Unterschied zwischen den beiden Newtonschen For- 
meln tritt auch in ihrer Verwendbarkeit fur praktiache Inter- 
polationszwecke deutlich hervor. Handelt es sich um die Inter- 
polation in der vollen Umgebung einer Stelle, so wird man zu 
der jiingeren Formel greifen. Kommt dagegen nur die halbe 
Umgebung in Frage, also das Gebiet rechts oder links von einer 
Stelle, 80 ist die altère Formel vorzuziehen. 

Anmerkungen. 

1. Zu S. 2. Man braucht in Wirklichkeit nur die ers t en 
Quotienten jeder Eeihe, d. h. a, |, A, l, und die erste Ordinate a, 

2. Zu S. 5. Das Bildungsgesetz der Faktoren von k, l, m ... 
tritt deutlicher hervor, wenn man die Formel so schreibt: 

k + xl+-m + — n+-^-^^o + j^^ 'p 

xHx'"'- \){x^ - 4) ^ , 

In dieser Fassung steht die Formel bei Lagrange (Astro- 
nomisches Jahrbuch von 1783). Stirling bringt durch Ver- 
einigung von je zwei Gliedern folgende einfachere Formulierung 
zustande (Methodus differentialis Newtoniana illustrata, 1719): 

, , 2lx + mx^ , An X •\- o x'^ a:* — 1 
^ + rT2 + 172 3Ti- 

^ px -{■ qx^ «*— 1 a;* — 4 
"^ ÏTi 3T4 5".^ T • . . 

3. Zu S. 6. Die Struktur der Formel tritt bei folgender 
Schreibung besser zutage: 

* + ^'+ 1.2 ''^^ 1.2.3 "+ 1.2.3-4 " 



X 



("-i)(--|)\(-T)(>-|)(--f) 



"^ 1.2.3.4.5 ^ "^ 1.2.3.4.5.6 ^ 



X 



(-i)(--T)(--f) 



^ 1.2.3.4.5.(5.7 ^ 
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Auch hier yereînigt Stirling je zwei Glieder und gelangt 
dadurch zu folgender Fassung: 

, . , . Sm + nx 4a;* — 1 , bo + px 4ic' — 1 4a;* — 9 
fl;-r^«'-r 4 ^ 2.3 ^ 4'' 2.3 4-5 

Iq + rx 4a;* - 1 4a;* - 9 ^ 4 a;* - 25 
"^ 4« ' 2.3 * 4.5 ' 6-7 "^ 

4. Zu S* 6. Um die Newtonschen Betrachtungen in 
tibersichtlicher Form durchzuftihren, muB man vor allen Dingen 
eine zweckmâBige Symbolik zugninde legen. 

Ftir einen ersten Diflferenzenquotienten 

?/« + ! - Vn 



ist die Bezeichnung 

ûblich, fiir einen zweiten Differenzenquotienten 

die Bezeiclinung 

iiir einen dritten Differenzquotienten 

[yn + i Vn + rVn+s] - [y» ^n + t Vn + î] 

die Bezeichnung 

L2/n 2/n + 1 2^n + 2 2^n + 3 J 

usw. 

Die altère Newtonsche Formel fiir die Parabel nter Ord- 
nung durch n + 1 Punkte [x^, y^), . . ., [x^, yj lautet dann 

2/ = 2/0 + (aï - ^o)[yoyi'] + (^ - aîo)(^ - ^i)[2/o2/iy2] + • • • 
. . . + (a; - a;,)(a; _ aj^) . . . (a; - a;„_i)[2/o2/i % • • • 2/J. 

Roger Cotes gibt hierfûr in seiner Abhandlung einen Be- 
weis, der an Einfachheit und Klarheit nichts zu wunschen iibrig 
lâfit Es genûgt aber eigentlich schon die zweite Aufstellung 
Cotes' (S. 15). Die Newtonsche Interpolationsformel hat nâm- 
lich die Eigentûmlichkeit, ihren Charakter zu bewahren, wenn 
man am SchluB ein Glied oder mehrere absôhneidet. LàBt 
man in der Formel nter Ordnung das letzte Glied fort. 
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80 entsteht die FoTmel {n — l)ter Ordnung. Um dies ein-' 
zusehen, braucht man sich nur folgendes klar zu machen: Wenn 

^0 + ^1 (^ - ^o) + ^2 (^ - ^o) {X-X^)+ ... 

• • • + ^«(^ - ^o) {x-x^) .,. [x- x^_^) 
fur Xç,. x^, . . ., x^ZM .%> 2/i, ' -> Vn wir<î> ^^ ^^^^ 

Co + Cj (x - Xo) + c^ (x - a:o)(^ - ^i) + - • • 

. . . + C„_i (X - X,)[X - Xj) . . . (X - X„_2) 

•fur a;^, Xi, . . ., x„_i zu t/^, 2/r • • •> yn-n ^^"^ ®^®^ ^^^ ^^^®^- 
scheidende Glied 

c^[x - Xq){x - x^) . , , [x - x^_^) 
gerade fur Xq, x^, . . ., x„_^ verschwindet. WeiB man nun aus 
Cotes' zweiter Aufstellung, daB 

^n = [% 2/i • • • 2/„] 

ist, 80 folgt au8 derselben Aufstellung 

^n-i =[2/o2/i • • • 2/„_i]- 
und weiter 

bis zu 

^1 =[2/o2/i]. ^0 = 2/0- 

Cotes' zweite Aufstellung und die sie sttltzende erste Auf- 
stellung bieten also einen sebr bequemen Weg zur âlteren Inter- 
polationsformel Newtons. 

Bedenkt man, daB es durch w + 1 Punkte nur eine Parabel 
wter (oder niedrigerer Ordnung) gibt, so wird sofort klar, daB 
das Newtonsche Polynom 

% + {x- x^,)[2/o2/i] + ... + (x- x^),.,[x- ir„_i)[2/o2/i---2/n] 

sich nicht ândern kann, wenn man die 7^ + 1 Punkte anders 
numeriert. Insbesondere muB also der Koeffizient von a;" der- 
selbe bleiben, d. h. [y^^ 2/i • • • 2/J verhàlt sich gegenûber allen Ver- 
tauschungen der n + \ Punkte indiffèrent. Dièse wichtige 
Symmetrieeigenschaft der Newtonschen Differenzen- 
quotienten bietet beim Eechnen mit ihnen groBe Erleichte- 
rungen. Merkwiirdig, daB weder Newton noch Cotes aus- 
drlicklich darauf hinweist. Von Newton freilich ist allgemein 
bekannt, daB er sich in seinen Publikationen immer groBe 
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Zurlickhaltung auferlegte. So lâBt er uns auch ûber den Weg, 
der ihn von seiner âlteren zu der jûngeren Interpolationsformel 
der Methodus dififerentialis ftihrte, ganz im Dunkeln. Die Be- 
merkung am Scblusse der zweiten Âafstellang (S. 3): ,,Âus 
diesen Âufstellungen lâfit sich das Folgende leicht herleiten^' 
ist ein sehr unzureichender Wegweiser. 

Wir môchten glauben, daB Newton ungefâhr so argumen- 
tiert hat, wie wir es jetzt darlegen werden. 

Erster Fall. 

Parabel durch eine ungerade Zabi von Punkten. 

Es bandle sicb z. B. um sieben Punkte mit den Nummern 
0, 1, 2, 3, 4, 5, 6. Wir scbreiben die altère Newtonsche Inter- 
polationsformel auf, wobei aber die 7 Punkte in einer anderen 
Folge genommen werden, und zwar in dieser: 

3, 2, 4, 1, 5, 0, 6. 

Sie entstebt ans der urspriinglichen dadurch, daB man zuerst 
den mittelsten Punkt 3 heraushebt und die so entstehende Lùcke 
abwecbselnd nacb links und nacb rechts yergrôBert, indem man 2 
und 4, 1 und 5, und 6 fortnimmt. Die Newtonscbe Inter- 
polationsformel lautet nun fur die neue Punktfolge: 

2/ = 2/3 + (^ - ^3) [2/3 y 2] + (^ - ^3) (^ - ^2) [% 2/2 2/ J 

+ {X- Xg) [X - X^) (X ~ Xj [2/3 2/3 2/4 2/1] 

+ {x-' Xg) {x - X.J,) {x - x^) [x - x{} [2/3 2/2 2/4 2/1 2/5] 

+ (X'- X3) [x - x.^) {x - x^) {x - x^) {x - X5) [% 2/2 2/4 2/1 2/5 2/0] 

+ {x- x^){x-x^){x-xJ{X''X^)[X'-x^){x-x^)[y^y.^y^y^y^y^,y^,], 

Legt man statt 3, 2, 4, 1, 5, 0, 6 die auf ganz âhnliche Weise 
entstehende Folge 

3, 4, 2, 5, 1, 6, 

zugrunde, so nimmt die Newtonsche Interpolationsformel die 
Gestalt an: 

?J = 2/3 + (^ - ^3) [2/3 2/4] + (^ - ^3) (^ - ^4^ [2/3 2/4 2^2] 

+ {x- X3) {x - x^) [X - X2) [y s ^^ y 2 2/5] 

+ (x - 0^3) {x - xj {x - x^) {x - x^) [y^ y^ % 2/5 2/i] 

+ {x- x^) [x - x^ {x - x^) [x - Xg) {x-x^) [2/3 2/4 y 2 2/5 2/1 2/0] 

+ (rc - aîg) (x - xj (a: - iCa) (a; - X5) (a; - X J (ir - Xe) [2/3 2/4 2/2 2/5 2/1 2^6 2/0] • 
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Aus beiden Pormeln ergibt sich durch Mittelbildung (Addition 
und Division durch 2) die jiingere Newtonsche Interpola- 
tionsformel fiir 7 Punkte. Sie lautet, wenn man die Sym- 
metrieeigenschaft der Newtonschen Differenzenquotienten be- 
rticksichtigt, folgendermaBen: 

+ (a: - x,){x ^ x,){x - X,) ^^-^-^-^-^ l [y^y^y^ y^l 

+ {x- x^)[x - x.^)[x - xj [x - ~^j [2/12/22/32/42/5] 

+ {x- X,) (X - X,) [x-x,] {X - X,) {x-^x,) \yoy.y.ysm y,\^y.y,ii.y^ 

4- (aï - x^) {x - x^) {x - xj (x - x^) [x -x-^)\x- -5l±^j [^^ y^ 2/2 2/3 2/^ 2/5 2/0] • 

Stellt man unter Zugrundelegung der ursprunglichen Folge 
0, 1, 2, 3, 4, 5, 6 das Schéma der Differenzenquotienten auf, 
80 erhâlt man: 

2/0 2/1 y 2 2/3 2/4 2/5 2/.; 

[2/02/1] [2/12/2] [2/22/3] [2/3 2/J [2/42/5] [2/52/0] 

[2/0-2/2] [2/1-2/3] [2/2-2/4] [2/3-2/5] [2/4 -2/0] 

[2/0 --2/3] [î/i'-î/J [2/2- -2/5] [2/3- -2/6] 

[2/0 •• • 2/4] [2/1 •• • 2/5] [2/2 •• • 2/0] 

[2/0----2/5] [2/i----2/g] 

[2/0 2/0] 

Man sieht dann deutlich, wie in der jiingeren Newton- 
schen Interpolationsformel die Mittellinie dièses Schémas aus- 
gezeichnet ist. 

Zweiter Fall. 

Parabel durch eine gerade Zahl von Punkten. 

Es handle sich z. B. um 6 Punkte mit den Nummern 0, 1, 
2, 3, 4, 5. Bevor wir die altère Newtonsche Formel anwenden, 
stellen wir âhnlich wie béim ersten Fall eine neue Reihenfolge 
her. Wir nehmen die beiden mittelsten Punkte heraus und er- 
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weitern die so entstehende Liicke abwechselnd oach links und 
recbts bzw. nach rechts oder links. Dadurch erhalten wir 

2, 3, 1, 4, 0, 5 bzw. 3, 2, 4, 1, 5, 0. 

Die zugehôrigen Newton schen Formeln 

y = 2/2 + [^- ^t) [î/2 J/a] + (^ - ^'ù (^ - ^3) [2/2 2/3 2/1] 
+ (rr - a;.) (x - x.^) [x - x^) [y., y^ î/i 2/ J 
-f (a; - a^j) (aï - aîg) (x-x{}(x- x^j [j/2 t/g 2/, î/4 î/u] 
+ (x - 3^2) (a;- aîg) (a; - aîj (a; - x^) (x - a;^) [y^ y^ y^ y^ y^ j/^] 

und 

2/ = 2/3 + (« - ^3) [2/3 2/2] + (^ - ^3)(^ - ^2)[2/3 î/2 Î/J 

+ (a: - aîg) (a; - x^) {x - x^) [1/3 7/^ y^ y{] 

-i'{x-x^)[x- X2) (x - x^) (x - x^) [î/g j/o î/4 2/1 2/5] 

+ (a- - ajg) (a: - x^) [x - x^] {x - ajj (a; - aîg) [j/g y^ 2/4 ^i ^s 2/o] 

werden nun durcb Mittelbildung zu einer neuen vereinigt, die 
eben die jûngere Newtonsche Interpolationsformel ftir 
6 Punkte ist. Sie lautet, wenn man wieder die Symmetrie- 
eigenschaft der Newtonschen Differenzenquotienten ausnutzt, 

+ {X- X-i) [x - X3) [x - -ÎLi-3.j [y^ y^ y^ j/ J 

+ {x- x,)(x - x,)(x - x,)ix - a;ji^l'Ay3Î^.l+Jî^,y_^»*.2/.l 

+ [x- x^) (x - Xj,) (a; - a;,) (a; - xj [x - ~—j [i/o Vi V-i y^ V^, y{\ ■ 
Das Schéma der Differenzenquotienten ist hier folgendes: 

I/o Vi Vi Va y* ys 

[yoz/J [i/i^a] [1/2 ys] [yjy*] [ytysJ 

[yo-yi] [yi-Vs] ^^yi-i/i] [ys-ys] 

[yo--ys] [yi-yJ [ya-s^s] 

[yo •••?/*} [.'/i-'-ysl 

[yo-'-'yi] 

Man sieht wieder deutlich, in welcher Weise die Mittellinie des 
Schémas bevorzugt ist. 
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Ftir den Fall gleicher Intervalle hat Stirling in der oben 
zitierten Methodus diflferentialis Newtoniana illustrata (1719) die 
hier dargelegte Dreiteilung des Newtonschen Verfahrens mit 
YoUer Klarheit entwickelt, ebenso im zweiten Teile seines 1730 
erschienenen Werkes Methodus differentialis sive Tractatus de 
summatione et interpolatione serierum infinitarum. Ftlr un- 
gleiche Teilintervalle beschrânkt er sich auf die altère Newton- 
sche Formel. Ûbrigens bespricht er auBerdem eine andere 
Interpolationsart, die sogenannte hyperbolische Interpolation, die 
tir sehr fernliegende Glieder geeignet ist, wo die Newtonschen 
Parabeln yersagen. 

n. Cotes. 

Roger Cotes (1652— 1716), ein Pfarrerssohn, zeigte schon 
in fruher Jugend groBes mathematisches Talent und erhielt in 
dessen Ausbildang jede nur môgliche Fôrderung, 1706 wurde 
er Professor der Astronomie und Physik in Cambridge. Newton 
schâtzte ihn sehr hoch und ubertrug ihm die Besorgung der 
2. Auf lage der Principia. AUgemein bekannt ist Cotes durch 
sein Theorem ttber regulàre Polygone, welches sein Vetter 
Robert Smith ans den Papieren des Verstorbenen gerettet hat. 
Smith gab 1722 die gesammelten Schriften Cotes' heraus unter 
dem Titel: Harmonia mensurarum sive Analysis et Synthesis per 
rationum et angulorum mensuras promotae. Accedunt alia 
Opuscula mathematica per Rogerum Cotesium. 

Unter den Opuscula oder Opéra miscellanea, wie der Sonder- 
titel auch heiBt, steht an zweiter Stelle die hier zum ersten 
Maie ins Deutsche iibertragene Abhandlung: De methodo differen- 
tiali Newtoniana. Cotes leitet darin die altère Newton sche 
Interpolationsformel ab, um sie auf das Quadraturproblem an- 
zuwenden. Er hat seine Resultate im AnschluB an die sehr 
knappen Newtonschen Bemerkungen in dem oben mitgeteilten 
Hilfssatz V des dritten Bûches der Principia gefunden (1707) 
und lernte Newtons Methodus dififerentialis erst 1711 kennen. 
Merkwtirdig ist es, daB Cotes den Unterschied zwischen seiner 
und Newtons Abhandlung hauptsâchlich darin erblickt, daB er 
die Koeffizienten der interpolierenden Funktion auf andere Weise 
bestimmt. Warum weist er nicht darauf hin, daB Newtons 
Methodus differentialis eine ganz neue Interpolationsformel bringt? 
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Cotes hat iibrigeDS auch iiber die jtingere Newtonsche 
Formel eine besondere Abhandlung verfaBt, die in den Opus- 
cula an dritter Stelle zu finden ist un ter dem Titel: Canono- 
technia sive Constructio tabularum per differentias. Lagrange 
macht in seiner Abhandlung iiber Interpolation, die 1778 in 
der Berliner Akademie gelesen wurde und 1783 deutsch im 
astronomischen Jahrbuch erschien, auf dièse Cotessche Ab- 
handlung und ihre Beziehung zu Newtons Methodus differen- 
tialis aufmerksam. Bei Moritz Cantor wird sie, wie mir 
scheint, gar nicht erwàhnt. Hâtte Cantor sie durchblâttert, so 
wàre er auf den Unterschied zwischen den beiden Newton- 
schen Formeln gekommen und batte sein ungenaues Référât 
ûber die Methodus differentialis korrigieren kônnen. Die Canono- 
technia ist ûbrigens etwas weitschweifig geschrieben und neben 
Newtons Methodus differentialis nicht von so groBer Bedeutung, 
daB sie in yorliegender Sammlung einen Platz verdiente. Stir- 
ling, dieser eifrige Interpret der Newtonschen Interpolations- 
methoden, der ûbrigens auch iiber den Newtonschen Katalog 
der Kurven dritter Ordnung eine erlâutemde Abhandlung schrieb, 
hat Cotes' Arbeiten iiber Interpolation nirgends zitiert. Auch 
in seinem Tractatus von 1730 spricht er nur ganz unbestimmt 
von „mehreren beriihmten Geometem, die nach Newton von 
der Beschreibung einer parabolischen Kurve durch irgendeine 
Anzahl von Punkten gehandelt haben*'. 

Anmerknngen. 

1. Zu S. 12. Vgl. den in der Eïnleitung der Anmerkungen zu 
Newton mitgeteilten Hilfssatz aus den Principia. 

2. Zu S. 13. Dièse Briefstelle (aus Newtons beriihmtem 
Brief an Leibniz bzw. Oldenburg vom 24.0ktober 1676) wird 
im Tractatus von Stirling (1730) sogar wôrtlich angefiihrt, wo- 
raus man sieht, wie groBes Gewicht die Anhânger Newtons 
auf dièse Datierung legten. Sie lautet in deutscher Uber- 
setzung: 

„. . . Wo die einfachen Reihen nicht leicht genug zu be- 
handeln sind, habe ich eine andere, nocl^ nicht verôffentlichte 
Méthode, durch die man sich dem Verlangten beliebig nâhem 
kann. Ihre Grundlage ist eine bequeme, leichte und allgemeine 
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Lôsung des folgenden Problems: Eine geometrische Kurve 
zu beschreiben, die durch beliebig yiele gegebene 
Punkte hindurchgeht. Euklides hat gelehrt, wie man 
einen Kreis durch drei gegebene Punkte beschreibt; es kann 
auch ein Kegelschnitt durch fiinf gegebene Punkte beschrieben 
werden und eine Kurve von drei Dimensionen durch sieben ge- 
gebene Punkte, 80 da6 ich die Beschreibung aller Kurven dieser 
Ordnung beherrsche, die nur durch sieben Punkte bestimmt 
werden. Dies lâBt sich sofort geometrisch ohne jede Zwischen- 
rechnung machen. Aber das oben genannte Problem ist von 
anderer Art. Und obwohl es beim ersten Anblick unangreifbar 
erscheint, so verhâlt sich die Sache doch anders. Es gehôrt 
nâmlich zu dem Schônsten, viras ich zu lôsen virûnschen 
môchte". 

Newton spricht hier von einer Kurve dritter Ordnung durch 
sieben Punkte. Er meint neun Punkte. Stirling ftigt (a. a. 0. 
S. 97), nachdem er die obige Briefstelle mitgeteilt hat, noch 
folgendes zur Erlâuterung hinzu: 

„ Newton lehrt im 58. Satze seiner Arithmefica universalis 
die Beschreibung einer konischen Parabel durch vier Punkte 
oder er lehrt vielmehr eine Méthode zur Auffindung der Glei- 
chung der Parabel, die durch vier gegebene Punkte hindurch- 
geht. Und nach derselben Méthode làBt sich eine Linie dritter 
Ordnung durch neun Punkte, eine Linie vierter Ordnung durch 
vierzehn Punkte beschreiben usw." 

^an sieht, wie taktvoU Stirling hier Newton korrigiert. 

3. Zu S, 16. Statt „dividierte Differenzen" oder Differenzen- 
quotienten sagt Cotes (wie Newton) vielfach kurz „Differenzen". 

4. ZuS. 17. Dies ist die altère New tonsche Interpolations- 
formel. In moderner Schreibweise lautet sie: 

2/ = % + (^ - ^o)[2/o2/i] + {x- x^){x - x^)[y^y^ 2/3] + • • • 

Die Koeffizienten sind aus dem Schéma der Differenzen- 
quotienten zu entnehmen. Aus der Newtonschen ergibt sich 
auch die beruhmte Taylorsche Formel durch einen Grenz- 
ubergang, den schon Taylor selbst in seiner Methodus incre- 
mentorum, wenn auch nicht bei beliebiger Intervallteilung, durch- 
gefuhrt hat (1715). 
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5. Zu S. 21, Das Beispiel der Binomialreihe entstammt 
dem ureigensten Gedankenkreise Newtons, der die Ausdehnung 
der Binomialformel auf beliebige rationale Exponenten durch 
Interpolation fand, Den Begriff der Interpolation verdankt 
Newton vermutlich Wallis, in dessen Arithmetica infinitorum 
(1655) davon die Rede ist. In Stirlings Methodus differen- 
tialis illustrata (1719), die wir bereits mehrfacb zitiert baben, 
steht dièses Beispiel auch. Stirling bringt noch mehrere an- 
dere schône Beispiele. Wir woUen hier wenigstens eins davon 
in modemisierter Fassung darlegen. 

Um 

u 

zu finden, interpoliert Stirling die Reihe 

X X 



X 



und findet fur y_^ oder arctgrc die bekannte Formel 

als deren Erfinder er Jacob Gregory nennt, der sie am An- 
fang des Jahres 1671 an CoUins mitgeteilt habe, von dem sie 
durch Oldenburgs Vermittelung an Leibniz gelangt sei („ad 
Leibnitium delata est^*). Man sieht, wie damais jede Gelegen- 
heit benutzt wurde, um auf den groBen deutschen Mathematiker 
zu sticheln. 

6. Zu S. 23. Dièse geometrische Quadratur wird erst klar, 
wenn man sie analytisch verdolmetscht. KS nimmt man am 
besten als Làngeneinheit. Die in stetiger Proportion stehenden 
(d. h. eine geometrische Reihe bildenden) Strecken KP, KS^ KT, 
KQ werden dann durch 

x — x^, 1, [x — x^)^\ {x — x^)~^ 
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gemessen, und es ist nach Cotes' Ângaben 

ù 

ki = i^^iTJ, ko = -^"•>' STO. 

o 4 

Der gesuchte Flâcheninhalt driickt sich, wie Cotes sagt, 
80 ans: 

AKPM=^{x - x^)KA + ^^"^^'^^ KE 

Was die Ordinaten KA, KE, KJ, KO anbetrifft, so ergeben 
sie sich, wenn man die vorliegende Kurve durch 

y=o, + c,{x-x,)+c,{x- x,Y + C3 (x - x,f 

approximiert C^,J C^ Cg, C3 sind dann gleich KA, KE, KJ, KO, 
Man erkennt nàmlich, daB sie die Werte von 



y - Vi ,, _ * - 



u = ^' , V = flir X = X, sind, 

X — Xi X — Xi ^ ' 

d. h. 



Hierzu vergleiche man die vierte Aufstellung Cotes' (S. 18). 

7. Zu S. 24. Gemeint ist die Methodus differentialis. 

8. Zu S. 25. Seine Quadraturentabelle bat Cotes an- 
scheinend erst im AnschluB an die Methodus differentialis von 
Newton berechnet. Dieser batte nur die zweite Formel an- 
gegeben. Cotes spricht sich nicht dartiber ans, wie er seine 
Eechnungen angestellt bat. Doch finden wir dies keineswegs 
auffàllig wie G au s s (vgl. dessen S. 94 mitgeteilte Selbstanzeige). 
Denn es ist wohl klar, daB Cotes einfach die Newton sche 
Interpolationsformel aufschrieb und dièse integrierte, ein aller- 
dings sehr umstândliches Verfabren. Bequemer geht es mit 
der Interpolationsformel von Lagrange, die Gauss benutzt. 
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m. Oauss. 

Gauss hat in der Théorie der genâherten Quadrataren 
dadurch einen wesentlichen Fortschritt erzielt, daB er sich von 
der Voraussetzung gleicher Intervalle frei machte. Dabei ergab 
sich die Frage, wie man in zweckmàBigster Weise die 
Teilung des Grundintervalls in nngleiche Stûcke zu wâhlen hat. 

Den besten Ûberblick tibèr den Inhalt der Gaussschen 
Arbeit und ihr Verhâltnis zu denen seiner Vorgànger gewinnt 
man aus der Selbstanzeige, die Gauss in den Gôttingischen 
gelehrten Anzeigen am 26. September 1814 erscheinen lieB. Sie 
sei hier voUstàndig abgedruckt: 

„Der kônigl. Sozietàt wurde am 16. September von dem 
Prof. Gauss eine Vorlesung eingereicht liberschrieben: 

Methodus nova integralium valores per approximationem 
inveniendi". 

Unter den verschiedenen Methoden zur genâherten Bestim- 
mung der Intégrale, oder wie es in der Sprache der âlteren 
Analysten hieB, zur genâherten Quadratur krummliniger Figuren, 
ist die Newton-Cotesische, welche sich auf die Interpolations- 
méthode grundet, eine der brauchbarsten. Newton batte eine 
Auflôsung der Aufgabe gegeben, durch eine beliebige Anzahl 
gegebener Punkte eine parabolische Kurve zu ziehen, deren 
immer leicht ausfiihrbare Quadratur dann nâherungsweise die 
Stelle der Quadratur der eigentlich vorgegebenen durch jene 
Punkte gehenden Kurve vertreten kann, und zwar desto genauer, 
je mehr Punkte man in Anwendung bringt. Newton batte es 
indessen bei dieser allgemeinen Andeutung bewenden lassen und 
nur gleichsam beispielsweise fiir den Fall von vier in gleichen 
Zwischenrâumen liegenden Ordinaten A, B, C, D den genâherten 
Flâchenraum zwischen der ersten und letzten, wenn deren Ent- 
fernung = R ist, durch 



(y^ + y^+I^+I^)^ 



angefûhrt. Cotes, welcher fiir sich, und noch ehe Newtons 
Schrift Methodus differentialis erschienen war, schon im Jabre 
1707 àhnliche Untersuchungen angestellt batte, wurde durch die 
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zierliche Form, in welcher Newton das Endresultat in obigem 
Beispiele dargestellt batte (pulcberrima'et utilissima regala nennt 
68 Cotes), bewogen, dièse Vorscbriften weiter und bis auf den 
Fall von 11 Ordinaten auszudehnen. Immer erscheint so der 
verlangte Flâchenraum in der Gestalt des Produkts der Basis, 
oder der Entfemung der âuBersten Ordinaten, in die Summe 
der durch bestimmte Zahlkoeffizienten multiplizierten Ordinaten, 
und zwar haben zwei gleich weit vom Anfang und Ende ab- 
liegende Ordinaten allemal gleiche Koeffizienten. Dièse Quadra- 
turkoeffizienten bis zu dem Fall von 11 Ordinaten gibt Cotes 
am SchluB der Abhandlung de methodo differentiali, welche 
einen Teil der Harmonia mensurarum ausmacht, ohne sich liber 
das Verfahren, wodurch er sie berechnet bat, weiter zu er- 
klàren. Vielleicht bat man es dieser anspruchslosen Kûrze, 
womit bloB das Endresultat dargestellt ist, zuzuscbreiben, daB 
dièse schône und zweckmâBige Méthode von den Analysten 
weniger gekannt und benutzt zu sein scbeint, als sie es verdient. 

Bei dieser Méthode liegt durchaus die Voraussetzung gleicher 
Abstânde zwischen den Ordinaten zu Grande. AUerdings scbeint 
beim ersten Anblick dièse Voraussetzung am einfacbsten und 
naturlichsten zu sein, und es war nocb nicbt in Frage gekom- 
men, ob es nicht dem ungeachtet noch vorteilbafter sein kônne, 
Ordinaten in ungleicben Abstânden zum Grunde zu legen. Um 
dièse Frage zu entscheiden, muBte zuerst die Théorie der 
Quadraturkoeffizienten in unbeschrânkter Allgemeinheit ent- 
wickelt und der Grad der Genauigkeit des Eesultats bestimmt 
werden. Es zeigte sich, daB die Bedingungen, wovon dieser 
Grad der Genauigkeit abhângt, von der Art sind, daB man die- 
selbe durch zweckmàBig gewâhlte Ordinaten in ungleicben Ab- 
stânden allerdings verdoppeln kann, so daB man mit einer be- 
liebigen Anzahl gehôrig gewâhlter Ordinaten ebenso weit reicht, 
als mit der doppelten Anzahl von Ordinaten in gleichen Ab- 
stânden. Dièse Untersuchungen, nebst der voUstândigen Théorie 
der zweckmâBigsten Auswahl der Ordinaten, der dabei anzu- 
wendenden Quadraturkoeffizienten und der Bestimmung des 
Grades der Genauigkeit, welchen dièses Verfahren gewâhrt, 
machen den Hauptinhalt der vorliegenden Abhandlung aus. 

Ans der kurzen Entwickelung der Théorie der Cotesischen 
Quadraturkoeffizienten, welche der Verf, vorausschicken zù 
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miissen glaubte, beruhren wir hier nur dasjenige, was den Grad 
der Genauigkeit betrififfc, welchen die dadurch gefundenen ge- 
nâherten Intégrale haben. Vor allem muB hier bemerkt werden, 
daB die Anwendbarkeit dieser Méthode, ebenso wie das Inter- 
polieren, auf der Voraussetzung beruhe, daB die Ordinaten 
innerhalb des zu quadrierenden Baumes sich durch eine konver- 
gierende Beihe darstellen lassen. Es sei x die Abszisse, y die 

Ordinate und das Intégral \ y dx werde von a; = ^ bis x — h 

verlangt Man fiihre statt x eine andere verânderliche GrôBe 
ein, indein man etwa 

X = g -\- [h — g)t 
oder auch 

^ ^ ^ij3 + h) + ^{h - g)u 

setzt. Hier muB also y sich durch Reihen wie 

a + at^- u' t'^ + a" t^ + usw. 
oder 

^ + ^' w + /S" w^ + r u^ + usw. 

darstellen lassen, die konvergieren, jene, wenigstens so lange <, 
dièse, so lange u nicht grôfier wird als 1. Man mag daher der 
Kûrze wegen den Koeffizienten a und /9' die Ordnung 1, den 
Koeffizienten ce" und /3" die Ordnung 2 usw. beilegen. Dies 
vorausgesetzt, wird gezeigt, daB die Fehler, denen man sich 
bei der Cotesischen Méthode aussetzt, zwar immer von einer 
hôheren Ordnung werden, je grôBer die Anzahl der zum Grunde 
gelegten Werte von y ist, jedoch so, daB eine ungerade Anzahl 
und die zunàchst grôBere gerade Anzahl immer Fehler von 
einerlei Ordnung hervorbringen. So ist fur drei Ordinaten der 
Fehler sehr nahe 

fiir vier Ordinaten nahe 



iîodann fiir -fiinf Ordinaten nahe 

1 



2688 



[h -g) à 



YI 
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und fiir sechs Ordinaten nahe 

= ^^('^-^)«" "«f- 

Man sieht bieraus, daB es im allgemeinen vorteilhaft sein wird, 
bei Ânwendung der Gotesiscben Metbode eine ungerade Anzabl 
Yon Ordinaten zu benutzen. 

Der Verf. gebt bierauf zu der allgemeinen Untersucbung 
liber, wo die Einscbrânkung^ dafi die Ordinaten gleicbe Abstânde 
Yoneinander baben, wegfâllt. Sind bier A^ A y Al' usw. die Werte 
Yon y, die entsprecbenden Werte Yon % bingegen a, a', a" usw. 
oder è, y, h" usw. die entsprecbenden Werte Yon u und ibre 
Anzabl n + 1, so wird das genâberte Intégral wiederum die 
Gestalt baben 

\}i-^g)[ItA + :R A -^H' A' + usw.), 

wo -R, E\ R", usw. Zablkoeffizienten sind, die unabbângig Yon 
der Funktion y bloB durcb a, a, a ' usw. oder durcb b, b\ h" usw. 
bestimmt werden. Die Untersucbungen des Verf. geben fûr 
dièse Bestimmung folgendes Résultat. Es sei 

T=[t-a)[t-a){t-a"). .. 

Aus der Multiplikation dieser ganzen Funktion Yon t, welcbe 
auf die Ordnung w + 1 steigt, in die unendlicbe Beibe, welcbe 

den Logaritbmen Yon — — - Yorstellt, nâmlicb 

^-1 + 1^-2 + 1^-3+ USW. 

ergebe sicb das Produkt T + T'\ so dafi T' die darin ent- 
baltene ganze Funktion Yon t bezeicbnet, so wie T" die ûbrige 
mit negatiYcn Potenzen Yon t ins Unendlicbe fortlaufende Beibe. 
Dies Yorausgesetzt, ergeben sicb die Quadraturkoeffizienten R^ 

R, R' usw., wenn man in -jm- fur t der Beibe nacb die Werte 

a, a\ a" usw. substituiert. Auf eine âbnlicbe und nocb etwas 
bequemere Art leitet man jene Koeffizienten aus ô, h\ h" usw. 
ab, indem man die Funktion 

i^u ^ h){u ^ V){u - V). . . 

durcb U, ihr Produkt in die unendlicbe Beibe 

M~^ + - w""^ + — w"^ -h iisw. 

3 

Kowalewski, Newton usw. 7 
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durch U' + U" bezeichnet (8o da6 U' die darin enthaltene ganze 
Funktion von u vorstellt), und dann in —ffr- fur u der Beihe 
nach die Werte 6, 6', V usw. substituiert Statt der gebrochenen 
Punktionen . „ , ^ J^ lassen sich auch ganze Funktionen 

a 1 a u 

von t und u finden, welche die Stelle von jenen vertreten 
kônnen, und ftir deren Bestimmung der Verf. eine allgemeitie 
Méthode entwickelt. 

Der Grad der Genauigkeit der Integrationsformel hângt 
nun von der Beschaffenheit der Beihe T" oder U" ab. Im all- 
gemeinen ist der Fehler zwar von der Ordnung n + 1 ; allein 
wenn von den ersten GUedern jener Beihen einige ausfallen, so 
wird der Fehler von einer hôheren Ordnung, so da6 wenn T" 
erst mit der Potenz ^-™ oder U" mit der Potenz w-"* anfângt, 
der Fehler von der Ordnung n + m vrird. 

Hieraus ergab sich nun, da6 insofem die Werte a, a', a" usw. 
oder ô, Vf h" usw. willkiirlich gewâhlt werden kônnen, dièse sich 
so bestimmen lassen miissen, da6 die ersten w + 1 Glieder von 
T" oder U" wirklich ausfallen, wovon die Folge sein wird, da6 
der Fehler der Integrationsformel auf die Ordnung 2n + 2 
kommt. Die Untersuchung schreitet demnach zu der Bestim- 
mung derjenigen FunktioDen T und U, fur jeden Wert von n, 
fort, wodurch der angegebenen Bedingung Geniige geleistet wird. 
Der beschrânkte Baum erlaubt uns nicht, in das Ëinzelne der 
Untersuchung hier einzugehen: wir bemerken also hier nur, daB 
dièse Funktionen ein sehr einfaches Fortschreitungsgesetz be- 
folgen und in genauem Zusammenhange stehen mit der Ent- 
wicklung der Beihen 

J u 4 ■ 

W^ + -^u'^ + —u"^ + -—u~'^ + usw. 

o D 7 

in kontinuierliche Bruche, die der Verf. in einer frûheren 
Abhandlung (Disquisitiones générales circa seriem inlinitam 

1 + ^^ X + etc.) gegeben hat. — Offenbar gibt demnâchst die 

Auf lôsung der Gleichung T = oder f/ = die Werte von 
a, a\ a" usw. oder b, b\ h" usw., und die Werte der Quadratur- 
ix « fûzienten werden nach den allgemeinen Eegeln bestimmt, die 
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in diesem Falle noch besondere Vereinfachungen vertragen. 
Ûbrigens werden allerdings in den meisten Fâllen sowobl die 
Werte von a, a\ a" usw. als die Quadraturkoeffizienten Irra- 
tionalgrôBen. Dies ist indes an sich sehr gleichgûltig, sobald 
nur ihre numerischen Werte ein fur allemal mit einem an- 
gemessenen Grad von Genauigkeit berechnet sind. Ist dies der 
Fall, so wird die Anwendung dieser Méthode auf irgendeine 
Anzahl von Ordinaten wenig oder gar nicht mehr Mûhe machen, 
als die Anwendung der Cotesischen Méthode auf eine ebenso 
groBe Anzahl, da hingegen letztere auf eine doppelt so groBe 
Anzahl angewandt werden mûBte, um ungefâhr dieselbe Genauig- 
keit des Résultats zu geben, wie erstere. 

Um fiir die Anwendung dieser neuen Méthode nichts zu 
wiinschen ubrig zu lassen, bat der Verf. noch die numerischen 
Werte von a, a', a" usw,, sowie von jB,^R', R" usw., auf 16 Dezi- 
malen berechnet, mitgeteilt, zugleich mit den Briggischen 
Logarithmen der letzteren auf lODezimalen, ailes bis zu dem 
Fall von sieben Ordinaten. In diesem letzten Falle wird der 
Fehler der Integrationsformel nahe 

— \ /v^V 

176 679 360 ' 

woraus man abnehmen kann, daB in den meisten in der Aus- 
libung vorkommenden Fâllen schon eine geringere Anzahl zu- 
reichen wird. Um die Anwendung der Vorschriften und ihre 
verhàltnismâBige Schârfe noch mehr zu versinnlichen, ist als 

Beispiel die Berechnuug von / -p ^ ^^^ a; = 100 000 bis 
X = 200 000 beigefiigt, wo schon bei der Anwendung von vier 
Werten der, Fehler nur -t^t^t^ttt^t^ des Ganzen ist, und bei einer 

Ob 000 000 

grôBeren Anzahl sich in den unvermeidlichen Fehlern verliert, 
die selbst die grôBeren Logarithmentafeln noch ûbrig lassen." 



Der Astronom J. F. Encke, ein Schiller von Gauss, macht 
ùber die Gausssche Abhandlung folgende intéressante Mitteilung 
(vgl. das Berliner astron. Jahrbuch von 1862 oder die Gesam- 
melten Abhandlungen Enckes von H. Gravelius, LBd. S. 102): 

„Gauss batte wahrscheinlich grôBere Hoffnungen auf dièse 
Art der Intégration gebaut, als er spâter verwirklicht fand. Er 
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beschâftigte sich zu der Zeit, als ich unter ihm studierte, mit 
seiner Abhandlung, von der die einzelnen Bogen mir zur Ab- 
schrift mitgeteilt wurden, die ich noch als ein wertvolles An- 
denken bewahre. Die Berechnung der Pallasstôrungen, die 
ebenfalls in dieser Zeit mit groBem Eifer von ihm verfolgt 
wurde, mochte mit diesen Untersuchungen in enger Verbindung 
stehen. Spâter hat er, soviel mir bekannt, keinen Gebrauch 
davon gemacht, wovon der Grund nicht schwer nachzuweisen 
sein dûrfte, und ûberhaupt mochte bei periodischen, grôBere 
Zeitrâume umfassenden Funktionen eine Anwendung dieser Mé- 
thode hôchst selten vorkommeD. Nur bei der Berechnung des 
numerischen Wertes eines einzelnen ganz bestimmten Intégrais, 
dessen Ermittelung auf anderem Wege allzu beschwerlich ist, 
wird man eine zweckmàBige Anwendung machen kônnen, wie 
Gauss es bei seiner Abhandlung in bezug auf das Intégral 

f-j^^ von X = 100 000 bis x = 200 000 als Beispiel aus- 

gefOhrt hat." 

Bemerkenswert ist noch, was Encke ûber eine praktische 
Priifung der Leistungsfàhigkeit der Gaussischen Méthode im 
Vergleich zur Cotesischen berichtet (a. a. 0. S. 124): 

„In den ersten fûnf Jahren der neuen Berliner Sternwarte 
beobachtete mein damaliger Gehilfe, der jetzige Direktor der 
Breslauer Sternwarte, Herr Prof. Galle, dreimal tâglich den 
Barometerstand und Thermometerstand zu den Zeiten, welche 
um den neunten Teil des ganze Tagbogens spâter als SoDnen- 
aufgang und friiher als Sonuenuntergang fielen, und zurzeit des 
nahen Mittags. Er erreichte dadurch dieselbe Genauigkeit, als 
wenn er fûnfmal tâglich in gleichen Zeitintervallen die Beob- 
achtungen angestellt hàtte, wie es in der Tat auch die Erfahrung 
bei der Herleitung des mittleren Thermo- und Barometerstandes 
ans seinen BeobachtuDgen, verglichen mit denen einer lângeren 
Reihe von Jahren, bestâtigt hat." 

Es handelt sich bei den Galleschen Beobachtungen um 
diejenige Abszissenwahl im Tagesintervall (0,1), die dem dritten 
Gaussschen Polynom T, also 

entspricht (vgl. Artikel 22, Nr. III der Gaussschen Abhandlung). 
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4- ist als Nâherungswert fiir a = 0,1127016653792583 benutzt. 
Die Mittagszeit wird durch die Wurzel a = 0,5 markiert. 
Y àes Tages vor Sonnenuntergang ist ungefâhr der zur dritten 

Wurzel a" gehôrige Zeitpunkt. 

Auffallend erscheint es, dafi man dièse Méthode der un- 
gleichen Intervalle sonst nirgends benutzt hat. 



Anmerkungen. 

1. Zu S. 26. r ist das Polynoin, welches die Lagrange- 
sche Interpolationsformel liefert Gauss scheint Lagranges 
Autorschaft nicht zu kennen oder nicht anzuerkennen (vgl. unsere 
2. Anmerkung zu Jacobi). Die Lagrangesche Interpolations- 
formel steht in den Leçons élémentaires sur les mathématiques 
données à TÉcole normale en 1795 (auch deutsch von Nieder- 
miiller, Leipzig bei Teubner). Dièse Vorlesungen wnrden 1812 
auf Lagranges Wunsch im Journal de l'Ecole polytechnique 
verôffentlicht (vgl. die FuBnote Serrets in seiner Lagrange- 
Ausgabe, Bd.Vn, S. 183). 

2. Zu S. 36. Wieder die Lagrangesche Interpolations- 
formel ohne Angabe des Autors. 

3. Zu S. 38. Man setzt am einfachsten 

T={t^a)T,, 
Dann folgk durch Differentiation 

also 

4. Zu S. 54. Es handelt sich um die beruhmte Arbeit 
ûber die hypergeometrische Reihe 

Die zitierte Formel lautet 
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log— - = 



-i^ 



2-2 

'-1TF'' 



'-U- 



4*4 „ 
7-9 



1 — U8W. 

5. Zu S. 55. Die Ausdrûcke W sind bis auf konstante 
Faktoren die sogenannten Leg en dr esche Polynôme. Sie lassen 
sich folgendermaBen darstellen: 

^^ """^ du- 
wobei c^ 80 zu wâhlen ist, da6 der hôchste Koeffizient 1 wird. 
Man erha.lt dièse Legendreschen Polynôme als EoefÊzienten^ 

wenn man 

1 



nach Potenisen von t entwickelt. Vgl. z. B. Heine: Handbuch 
der Kugelfunktionen, wo tibrigens im 2. Bande die Gausssche 
Théorie voUstândig dargestellt wird. 

6. Zu S. 57. Die zitierte Formel lautet 

i^(«, /3 + 1, / +.1) - F{a, ^, y) = ^r^^^ F{u + 1,^ + 1,7 + 2). 

IV. Jacobi. 

Der Kerngedanke der aus Bd. 1 des Crelleschen Journals 
hier mitgeteilten deutschen Arbeit Jacobi s besteht darin, da6 die 
Gausssche Fragestellung auf folgendes Problem reduziert wird: 

Ein Polynom n-ten Grades zu finden, das mit jedem Polynom 
niedrigeren Grades ein Produkt liefert, dessen Intégral (lângs 
des Grundintervalles) verschwindet. 

Jacobi hait es fur unnôtig zu sagen, daB dièses Problem 
mehr als 40 Jahre vor ihm von Laplace und Legendre gelôst 
war, und zwar bildet das n-te Legendr esche Polynom 

. d-\{x- nr(x-br\ 
^ ' daf" 

die Losung, wenn a,, .b das Grundintervall ist 
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Hâtte Jacobi die Anknûpfung an dièse Vorarbeiten nicht 
yerschmâht, so wâre sein Beweis auf die Hâlfte des Umfanges 
zusammengeschnimpft. 

Wir wollen noch bemerken, daB man das Gausssche 
Problem aucb als reines Interpolationsproblem bebandelt bat. 
Die Âufgabe lautet dann: 

Wie muB die Newtonscbe Parabel (n— l)-ter Ordnung 

gew&blt werden, damit sie die gegebene Eurve 

y - m 

moglicbst genau im Intervall a,,, h darstellt? 

FaBt man die Genauigkeit im Sinne der Metbode der 
kleinsten Quadrate auf, so wird gefordert, daB das Intégral 



^[f[x)--P{x)Ydx 



ein Minimum sei. Nacb der âlteren Newtonscben Formel 
(vgl. Anmerkung 4 zu Newton) ist aber 

Das Intégral nimmt also folgende Gestalt an 
h 

a 

oder nacb dem Mittelwertsatz 

C2/oyi---2/n-i^?-/{(aî-aîo)-"(iC-a!n-i)P<^a;. 

a 

Uber den erstôn Faktor dièses Produktes lâBt sicb nicbts 
sagen, da er von der unbekannten Funktion f[o^ abbângt. Da- 
gegen ist der zweite Faktor ganz unabbângig von f[x). Es liegt 
nabe, durcb môglicbste Verkleinerung des zweiten Faktors eine 
Herabdrûckang des Produktes zu versucben. Man kommt bier- 
durcb auf die Frage: 

Wann ist das Intégral 
h 

J{[x-Xq)...[x- x^_,^)Ydx 

a 

ein Minimum? 
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Die Antwort lautet: x^, x^, ..., a;„_i mtissen die Wurzeln 
des n-ten Legendreschen Polynoms sein. Das deckt sich aber 
mit dem Gaussschen Ergebnis. 



Anmerkungen. 

1. Zu S. 69. Gemeint ist das inehrfach zitierte Lemma Y 
des 3. Bûches der Principia. 

2. Zu S. 70. Die Partialbruchzerlegung wurde 1702 von 
Leibniz gefunden. Durch Multiplikation mit dem Nenner q>{x) 
ergibt sich aus ihr die Lagrangesche Interpolationsfdrmel, die 
ûbrigens neuerdings auch Waring zugeschrieben wird. 

3. Zu S. 74. Der Zusammenhang zwischen den Integralen 

l q>{x)dx, l x(f{x)dx, . . . , l x^~'^(p(x)dx 
und dem w-fachen Intégral 



/' 



(p [x) d x^ 
lâBt sich am kûrzesten durch die Formel 

n X 



ausdrûcken. 

4. Zu S. 75. Man vermiBt hier den Hinweis auf Legendre. 
Jacobi hat ihn in einer spâteren Arbeit „Uber eine besondere 
Gattung algebraischer Funktionen, die aus der Entwicklung der 
Funktion (1 - 2a;aj + ^c^j-v. entstehen'' (Crelles Journal, Bd.2) 
nachgeholt. Er zitiert dort Legendres einschlâgige Abhand- 
lungen aus den Jahren 1784 und 1789 und entwickelt ver- 
schiedene neue Eigenschaften der Legendreschen Polynôme, 
so die Darstellbarkeit des wten Polynoms X^"*^ als wte Ablei- 
tung. Zum SchluB heiBt es dann: „tJbrigen8 ist die Funktion X^"* 
dieselbe, welche Gauss in seiner Abhandlung Methodus nova 
integralium valores etc. mit U bezeichnet. . . ." 

7 

«-2 
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